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Abstract 

In 1967 M.C. Gutzwiller succeeded to derive the semiclassical expression of the quantum 
energy density of systems exhibiting a chaotic Hamiltonian dynamics in the classical limit. 
The result is known as the Gutzwiller trace formula. 

The scope of this review is to present in a self-contained way recent developments in 
functional determinant theory allowing to revisit the Gutzwiller trace formula in the spirit 
of field theory. 

The field theoretic setup permits to work explicitly at every step of the derivation of 
the trace formula with invariant quantities of classical periodic orbits. R. Forman's theory 
of functional determinants of linear, non singular elliptic operators yields the expression 
of quantum quadratic fluctuations around classical periodic orbits directly in terms of the 
monodromy matrix of the periodic orbits. 

The phase factor associated to quadratic fluctuations, the Maslov phase, is shown to be 
specified by the Morse index for closed extremals, also known as Conley and Zehnder 
index. 

Key words: Semiclassical theories and applications (PACS:03.65.Sq), Classical and 
semiclassical techniques (PACS:11.15.Kc), Path-integral methods (PACS:31.15.Kb). 
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Introduction 



The Gutzwiller trace formula is the semiclassical expression of the energy density 
of a quantum system [81,82,83,84,87,88]. Despite of the mathematical difficulties 
with which it is intertwined, it represents the cornerstone of the present under- 
standing of the manifestations of the quasi- stochastic nature of the trajectories of a 
generic, chaotic, classical system in the properties of the quantum mechanics asso- 
ciated to it by Bohr's correspondence principle. 

The trace formula states that the energy spectrum of a non relativistic quantum sys- 
tem is given by a series over all the periodic orbit of the classical system. The dis- 
covery of the trace formula revived the intuition of Bohr and Sommerfeld [28,148] 
in the early days of quantum theory that the energy spectrum of a generic quantum 
system can be expressed in terms of the invariant properties of the periodic orbits 
of the corresponding classical system. The derivation of the trace formula is highly 
non trivial due to the singular nature of the classical limit of quantum mechanics. 
Feynman's path integrals offer the physically most transparent route to take the 
limit. Indeed path integrals were the starting point of Gutzwiller's analysis. At the 
same time, the very discovery in the late sixties of the trace formula corresponded 
also to a major progress in the understanding of real time path integration. Namely 
Gutzwiller realised that the known Van Vleck-Pauli-Morette [1 17,127,155] approx- 
imation of the semiclassical propagator had to be corrected by the introduction of 
a further phase factor in correspondence of conjugate points of the classical trajec- 
tory. The phase factor turned out to be specified by the Morse index for open ex- 
tremals [118]. The occurrence of the Morse index in the semiclassical propagator 
is experimentally observable. Namely, the phase of the propagator enters the Bohr- 
Sommerfeld quantisation condition and its generalisations as they emerge from the 
trace formula. 

The role and significance of the propagator phase were further clarified by the en- 
suing fundamental investigations of the semiclassical limit carried out by Maslov 
[109], Voros [158] and Miller [111]. 

The works of Gutzwiller and Maslov laid also the basis for the development of func- 
tional semiclassical methods in field theory. Gutzwiller's idea to compute atomic 
spectra without resorting to the construction of wave functions is of great advan- 
tage in field theory where state functionals are exceedingly complicated objects. 
Dashen, Hasslacher and Neveu [43,44,45] used the trace formula as a paradigm in 
their ground-breaking investigations of particle spectra of field theories. 
Since the second half of the seventies the kinds of semiclassical approximations 
underlying the trace formula became familiar tools of theoretical physics. Among 
the early applications one can mention the monopole quantisation in the Georgi- 
Glashow model by 't Hooft [152] and Polyakov [130] and the development the 
instanton formalism of Belavin, Polyakov, Schwartz and Tyupkin [17] and 't Hooft 
[153]. 
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The semiclassical approximation of path integrals is an infinite dimensional imple- 
mentation of stationary phase methods. Fluctuations around the field configuration 
dominating the path integral in the semiclassical limit are, within leading order, 
generically quadratic. The resulting infinite dimensional Fresnel integral brings 
about the exigency of computing the square root of the determinant of the self- 
adjoint operator governing the quadratic fluctuations. Thus, a natural mathematical 
counterpart of path integral semiclassical methods is the theory of functional deter- 
minants. 

Non-relativistic systems with strictly positive definite kinetic energy offer an intu- 
itive framework to understand the mathematical issues interwoven with the evalu- 
ation of functional determinants. Since the quantum dynamics can be described by 
means of configuration space path integrals, the self-adjoint fluctuation operators 
emerging in the semiclassical approximation are of Sturm-Liouville type with spec- 
trum bounded from below. The finite number of negative eigenvalues, the Morse 
index, by fixing the value of the phase of the quadratic path integral, provides a nat- 
ural definition of the winding number of the functional determinant in the complex 
plane. 

An intrinsic interpretation of the result is achieved by embedding the fluctuation op- 
erator in a family the elements of which are related by homotopy transformations. 
Thus, the phase of the functional determinant is specified by the winding number 
of the determinant bundle obtained by connecting the elements of the family to a 
positive definite operator with zero Morse index. Homotopy transformations may 
clip the explicit form of the operator as well as the boundary conditions obeyed by 
its functional domain of definition. A satisfactory theory should then be able to re- 
solve the dependence of functional determinants versus the parameters entering the 
operator and the boundary conditions. A powerful result was recently obtained by 
R. Forman [67] for functional determinants of a rather general class of elliptic oper- 
ators. In the context of non-relativistic quantum mechanics Forman's result allows a 
particularly compact and intuitive expression of the semiclassical approximation of 
path integrals with general boundary conditions. The result has the further merit to 
provide a topological characterisation of extremals of classical mechanical action 
functional independently of the adoption of a Lagrangian or Hamiltonian formal- 
ism. In consequence of Forman's result, it is possible to establish the equality of 
the functional determinants of the self-adjoint operators associated to the second 
variation around a classical trajectory in configuration and phase space. In this lat- 
ter case quadratic fluctuations are governed by Dirac type operators the spectrum 
whereof is unbounded from below and from above. The proof of index theorems for 
such Dirac operators is based on the construction of an infinite dimensional Morse 
theory ultimately relying on spectral (Fredholm) flow theorems (see [137,143] and 
references therein; the particular case of open extremals was previously studied in 
the physical literature in [101,103,114] ). In the periodic case, the resulting index 
is named in the mathematical literature after Conley and Zehnder who showed its 
role in the proof of existence statements for periodic solutions of time periodic and 
asymptotically linear Hamiltonians [35]. The equality of the functional determi- 
nants entails immediately that the Conley and Zehnder index coincides with the 
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Morse index for closed extremals whenever both of them are defined. The result is 
proven in the mathematical literature by different, less intuitive methods [3]. In the 
physical literature an equivalent result was obtained by Sugita in [150]. 

The modem theory of functional determinant provides considerable insight in the 
Gutzwiller trace formula. 

In the traditional derivation the resort to path integration is restricted to the con- 
struction of the quantum propagator. This is equivalent to applying functional de- 
terminant theory only to fluctuation operators with Dirichlet boundary conditions 
exploiting a classical result of Gel'fand and Yaglom [72]. The semiclassical trace 
operation is carried out in a separate step. The re-parametrisation invariance [130] 
of closed quantum paths contributing to the density of states is broken in the pro- 
cess of the stationary phase approximation. In consequence, classical periodic or- 
bits do not appear at intermediate steps. They are shown to fully characterise the 
final result only after nontrivial algebraic manipulations which require a good deal 
of a priori physical insight. More seriously the canonical invariance of the trace, al- 
though physically expected, needs to be proven apart. The proof relies on an highly 
non trivial mathematical apparatus developed by Amol'd [8] and investigated in the 
physical literature by Littlejohn, Creagh and Robbins [39,105,140]. 
The situation is much simpler if the stationary phase approximation is carried out 
directly on the loop space where trace path integrals have support. Periodic func- 
tional determinants completely specify the contribution of quadratic fluctuations. 
The topological invariance of the phase factors is guaranteed by the Morse index 
theory for closed extremals [118]. The topological formulation of Morse index the- 
ory of Bott [29] and Duistermaat [56,40] renders then a priori available an arbitrary 
number of equivalent prescriptions for the explicit evaluation of the index. 

The aim of the present report is to review at a level as elementary as possible the 
general mathematical methods needed for a pure path integral derivation of the 
Gutzwiller trace formula. The systematic application of these methods to the en- 
ergy density of a non relativistic quantum system in the semiclassical limit is it- 
self a novelty. The merits of the path integral formulation are particularly evident 
when the Lagrangian of the classical system is invariant under extra continuous 
symmetries beyond time translations. In the path integral formalism, such situa- 
tions present no conceptual difference from the case when the energy is the only 
conserved quantity. In all cases trace path integrals are handled in a canonically 
invariant way by means of the Faddeev-Popov method [60,61]. 

The report is organised as follows. In chapter 1 the definition of the quantum den- 
sity of states is recalled. The expression of the semiclassical approximation found 
by Gutzwiller is then qualitatively discussed. Chapter 2 focuses on functional de- 
terminant theory and Forman's identity. Since the motivation is the semiclassical 
approximation of path integrals, the probabilistic interpretation of these latter ones 
is recalled in the first section of the chapter. Further details together with an out- 
line of the recent rigorous proof [5] of the existence of the covariant representation 
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of the path integral measure are given in appendix C. Chapter 3 summarises ba- 
sic results of Morse index theory [118,112]. The goal is to motivate the use of 
the symplectic geometry and differential topology concepts applied in the modern 
formulation of the theory [56,40,143,4] and to emphasise their physical relevance. 
The theory is illustrated with examples. Finally, it is shown how to derive from the 
general formalism practical prescriptions to compute the phase factors intervening 
in the Gutzwiller trace formula. Chapter 4, the last, deals with implementation of 
the Faddeev-Popov method for trace path integrals. The Gutzwiller trace formula 
is then proven to follow immediately from the application of the method. 
The material is expounded in each chapter in an as far as possible self -consistent 
way in order to allow independent reading. The main text is complemented by ap- 
pendices summarising basic concepts of stochastic calculus and classical mechan- 
ics. 

The author's wish is that this may serve the reader interested not only in the Gutzwiller 
trace formula but also in understanding methods of path integration of general use 
in theoretical physics. 
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1 The density of states in quantum mechanics 



The aim of the present chapter is to set the scene for the path integral methods which 
will be illustrated in the rest of this work. Some basic concepts of non relativistic 
quantum dynamics are shortly recalled in order to derive the relation between the 
density of states and the propagator. The relation is the starting point for the semi- 
classical approximations which led Gutzwiller to write his trace formula. Finally, 
the qualitative features and the significance for quantum mechanics of the trace 
formula are shortly discussed. 



1.1 From the Schrddinger equation to the energy density 



Non relativistic quantum mechanics for a spin-less particle is governed by the 
Schrodinger equation 



i h OtiP — Hijj = 

(1) 

where is the wave function, H the Hamiltonian operator and h Planck's constant. 
This latter is a pure number if absolute units are adopted 



[P]^-[Q]^1 (2) 

Q and P being respectively the spatial position and the canonically conjugated mo- 
mentum variables. The wave function t/j is in general a complex valued function. It 
specifies the probability amplitude to observe the system in the position Q at time 
T: the modulus squared \ip\'^ defines the probability density of the same event (see 
for example [53,64,98,141]). 

In typical situations a quantum particle of mass m interacts with an external poten- 
tial U and with a vector potential Aa. Most often the physical configuration space 
Wl is modeled by the Euclidean space W^, eventually represented in non Cartesian 
coordinates. A slight generalisation is to imagine 971 to be a Riemann manifold 
which is either compact or and has time independent metric Qap-^^ such a case 
the Hamilton operator is 



ff^:=ij«<'(^i±iii(5l + i£)v, + aV. (3) 
2 m 

where Aa and U are functions of Q and eventually T and (?°'^ is the inverse of gag. 
In (3) the Pa's denote momentum operators in position representation. They act on 
all objects to their right as covariant derivatives 
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Pa := -ihVa 
and therefore they obey the Leibniz rule 



(4) 



P„A"V = {PaA")lP + A"P„^ (5) 

The explicit form of the covariant derivative is fixed by the compatibility condi- 
tion with the metric assigned to SDT. These geometric concepts are recalled in ap- 
pendix A. 

The Hamilton operator acts on the space Lgjj of square integrable functions on SDT 
such that the Schrodinger equation is self -adjoint with respect to the scalar product 



g := det{g,,0} (6) 

If the potentials and U are time independent it is possible to introduce the sta- 
tionary Schrodinger equation 



E^Pe-H^e = Q (7) 

describing a quantum phenomenon occurring at a constant energy E. For bounded 

physical systems, the stationary equation admits solutions only if the energy as- 
sumes discrete, quantised, values. The energy levels are in such a case labelled by 
a set of integers referred to as quantum numbers. 

The imaginary factor appearing in the Schrodinger equation (1) ensures invariance 
under time reversal of the dynamics of isolated quantum systems. Nevertheless it is 
useful to introduce the forward time evolution operator or propagator. The kernel 
K of the propagator is specified by the solution of 



[ihdT - H)K{Q,T\Q' ,T') = -i5{T -T')5{Q - Q') , T > T' 

(8) 

K{Q,T\Q',T')^0, T<r 

The propagator is the inverse of the Hamilton operator for Cauchy boundary con- 
ditions in time. The propagator owns its name because it governs the propagation 
forward in time of wave functions 



m,T)^j d''Q'4^)K{Q,T\Q',T')^{Q',T') (9) 

JOT * 

If the Hamiltonian is time independent and the energy spectrum discrete the prop- 
agator is amenable in the sense of to a series over the eigenfunctions of the 
stationary Hamilton equation 
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KiQ,T\Q',T')^^iT-r)Y.^P^{Q)rniQ')e-'^'^ (10) 

n 

where 6{T — T') is the Heaviside step function 

f 1 if T > T' 

e{T-T') = { ~ (11) 

[ if T < T' 

Analogous representations exist for continuous and mixed spectra. 

The energy spectrum of the quantum system can be computed from the Fourier- 
Laplace transform of the propagator. The latter is given by the analytic continuation 
in the upper complex energy-plane 

Jo fi ^ Z — t,n 

Z eC, ImZ > (12) 

The result is proportional to the kernel of the Green function G, the resolvent of the 
time independent problem 

{Z -H)G{Q\Q',Z) = 5{Q-Q'), ImZ>0 

poo rj'~r yrj, 

G{Q\Q',Z)^ e^-K{Q,T\Q\0) (13) 

Jo in 

The poles of G are located on the real axis and coincide with the energy levels of 
the quantum system. Using the Plemelj identity [162] 



= P.y.-^ -inS{E-E^) (14) 

E=RcZ — J^n 

where P.V. denotes the integral principal value, the announced relation between 
the energy density of the system and the trace of G is found 



lim 

ImZjO Z - En 



p{E) = - lim Im- / d'Q 49iQ)G{Q\Q, Z) 



(15) 

E=ReZ 



Physically (15) means that the energy density is a scalar quantity in non-relativistic 
Quantum Mechanics, independent of the representation of the Hilbert space. 

Finally, in terms of the propagator the energy density reads 
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Jo 



•0° dT 




d''Q^f^)K{Q,T\Q,0) (16) 



inh 



E=ReZ 



1.2 The semiclassical limit and the Gutzwiller trace formula 

The relation between the trace of the propagator and the energy density provides a 
way to extract information on the spectrum of a quantum system without solving 
the stationary Schrodinger equation (7). Indeed the propagator can be regarded as a 
more fundamental object than the Schrodinger equation. According to Feynman's 
formulation of quantum mechanics [62,64] the propagator is the "sum" 



extended over all paths in configuration space connecting Q to Q' in the time in- 
terval \T' , T]. Each path contributes with a phase specified by the action functional 
S. The precise meaning of the sum will be recalled later. Here it is interesting to 
observe that classical mechanics is recovered in the singular limit of h tending to 
zero. The limit is directly meaningful in absolute units. In general units it will cor- 
respond to the vanishing of an overall adimensional parameter obtained from the 
Planck constant times other invariant quantities of the physical system. 
Bohr's correspondence principle requires the recovery of classical mechanics when 
h tends to zero while all the other parameters are held fixed. The existence of the 
limit entails the possibility to investigate a range of phenomena for which quantum 
effects are weak by means of an asymptotic expansion around the classical limit, 
the semiclassical expansion. 

Broadly speaking the semiclassical approximation is expected to apply to the de- 
scription of phenomena occurring at energy scales large in comparison to the mean 
energy level spacing. Typical examples are encountered in the context of meso- 
scopic physics where the investigation of highly excited states of atoms as well as 
the transport properties of solid-state devices are amenable to semiclassical meth- 
ods. However, the domain of applicability extends even to the ground state of cer- 
tain classes of quantum systems. The rescaling of the action functional S often 
evinces that the small h limit is equivalent to a small coupling regime of some non- 
linear term [34]. This observation has proven to have far reaching consequences in 
the investigation of tunneling phenomena where analytic perturbation theory is not 
available, both in systems with a finite number (Quantum Mechanics) and infinite 
(Field Theory) number of degree of freedoms. Finally loop expansions around a 
"vacuum" state can be ordered in powers of h [165]. 

The main result of semiclassical methods is to express quantum observables in 
terms of classical objects. One of the striking differences between classical and 



K(Q,T\Q',r) 




(17) 



paths 
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quantum dynamics is that the latter gives a linear evolution law for the probability 
amplitudes. In Classical Mechanics a generic, non-linear, Hamiltonian exhibits a 
chaotic dynamics: exponential sensitivity to initial conditions in bounded regions 
leading to stretching and folding in classical phase space. In Quantum Mechanics 
there is sound evidence that the evolution of time dependent observables, broadly 
speaking quantities related to the squared absolute value of probability amplitudes, 
exhibits no chaotic behavior [20] . However, there are many experimental evidences 
(see [30] for review) supporting the conjecture that the energy spectrum may re- 
flect the properties of a classically chaotic motion. The inference is motivated by 
the correspondence principle. The quantum energy levels are determined by the 
existence of eigenfunctions of the Hamilton operator. The correspondence princi- 
ple leads to associating them to the invariant sets of the classical dynamics. For a 
generic autonomous system invariant sets are the energy surface, the tori produced 
by eventual extra symmetry of the theory and classical periodic orbits. 
The role of classical invariant sets in quantisation is understood in the particular 
case of classically integrable systems. In many cases integrability follows from the 
possibility to separate variables ^ . A classical Hamiltonian system is separable in 
d dimensions if, including the Hamiltonian function H, there are d independent in- 
tegrals of the motion Hi,..., Ha-i) Poisson commuting with each other. The 
set of first integrals (7Y, Ha-i) is then said to be in involution. Using these 

d first integrals one can introduce action-angles coordinates {Ai, Ad, 9i, ...,9d) 
which are canonical coordinates such that the Hamiltonian H. = 7i(Ai, Ad) and 
the other first integrals become functions of the action variables alone. If a phys- 
ical system is classically separable, semiclassical quantisation gives a simple rule 
for the spectrum of a complete set of observables. Essentially this is a generalisa- 
tion of the Bohr and Sommerfeld rule (see for example [28,148,98,30]) and it is 
known as the Einstein-Brillouin-Keller or EBK quantisation rule [58,92]. The EKB 
quantisation predicts the quantisation of the action variables 



with the Tij's and fi/s specifying the set of quantum numbers of the system [98]. 
The relevance of (18) is mainly conceptual. It establishes a remarkable connection 
between the occurrence of classical periodic orbits with quantum spectra. Namely, 
the quantised values of the action variables in (18) correspond to closed curves on 
the tori defined by the constants of the motion. From the practical point of view, 
the Schrodinger equation always separates when the classical problem is separable 
[87] rendering thus available the exact expression of the energy spectrum. 
In generic classical systems there are no constants of the motion other than the 
energy. Moreover it is known that in a chaotic system the number of primitive 
periodic orbits proliferates exponentially with the period T. The phenomenon has 

^ A famous counter example is the Toda lattice [89,65] (see also [80] for review). The 
quantisation of the Toda lattice is investigated in [85,86]. 




(18) 
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no correspondence in the observed densities of energy levels. 
A further difficulty, pointed out by M. Berry, towards the understanding of how 
classic chaotic behavior translates into quantum spectra arises from the infinite time 
limit which is intertwined with both concepts [20,22]. The very definition of typical 
indicators of chaotic behavior like of Ljapunov exponents, Renyi entropies (see for 
example [9,15,42]) and algorithmic complexity [66] requires infinite sequences of 
data. Such data are provided by the same infinite time limit which defines invariant 
states. However, the correspondence principle does not guarantee that the classical 
limit will be interchangeable with the infinite time limit and counterexamples are 
known [22]. 

The difficulties listed above stress the relevance of the breakthrough achieved in a 
series of impressive papers [81,82,83,84] by M. Gutzwiller with the discovery of 
the trace formula which now bears his name. Gutzwiller's trace formula yields a 
general semiclassical expression for the quantum energy density. In the case when 
the energy is the only conserved quantity, the trace formula takes the form 



p{E) ^J:5(e- ^ 

n 

±^^s{E-n{p,Q)) + im y: 1 (19) 



The semiclassical density of states is seen to comprise two terms of different na- 
ture. 

The first term is a microcanonical average of the classical Hamiltonian, associated 
by Bohr's correspondence principle to the Hamilton operator the spectrum whereof 
is sought. The existence of a similar contribution has been known for long time in 
atomic physics from the Thomas-Fermi approximation. The microcanonical aver- 
age brings in a smooth background dependence of the energy density on E. 
The second term, Gutzwiller's genuine achievement, consists of a formal series 
ranging over all classical primitive periodic orbits of finite period To and their repe- 
titions r. Each orbit is represented in the series by a complex function of the energy. 
The phase has an essential singularity in h proportional to the reduced action W of 
the periodic orbit [98]. The reduced action depends continuously on the energy and 
therefore brings in strong oscillations in the energy density. The phase also receives 
contribution from a topological invariant of classical orbits usually referred 

to as the Maslov index. It carries information about the structural stability of the 
dynamics linearised around the orbit. 

The amplitudes of the orbit contributions depend upon the monodromy matrix M 
of each primitive orbit. The monodromy matrix governs the linear stability in phase 

space of the periodic orbit over one period. It enters the trace formula through the 
inverse square root of the absolute value of the determinant of \2d — M restricted 
to the eigendirections transversal to the orbit. Thus, the contribution of the more 
unstable orbits is exponentially damped 
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^|det^(l2d-M(r))| 



exp < — } (20) 



the decay rate being specified by tlie Kolmogorov-Sinai entropy [147,9,15,125]: 
the sum of the positive Ljapunov exponents of the orbit. 

Periodic orbits are therefore seen from the Gutzwiller trace formula to affect the 
spectrum both individually, through the essential singularity in h, and collectively. 
The collective contribution gives rise to major physical and mathematical difficul- 
ties. For fixed values of the energy the number of periodic orbit is infinite. More- 
over, in a chaotic system the number of periodic orbits proliferates exponentially 
with the period T and growth rate given by the topological entropy hx'. 



# (periodic orbits) ~exp{/iTT}, T t oo (21) 

The topological entropy is the Renyi entropy of order zero [147,9,15,125]. If one 
assumes that on average the topological and Kolmogorv-Sinai entropies are equal, 
the diminishing amplitude of orbits of period T is dominated by their proliferation. 
Thus, the series consists effectively of terms with exponentially growing ampli- 
tudes. The estimate indicates that a literal interpretation of the Gutzwiller trace 
formula is problematic. Nevertheless experimental and numerical evidences exten- 
sively reviewed in [30] support the existence in some mathematical sense of a semi- 
classical approximation to the density of states related to the trace formula above. 
In particular, it has been conjectured that convergence may arise on the basis of 
the topological organisation which is often observed in the occurrence of periodic 
orbits. Longer orbits should be "shadowed" by shorter ones begetting in this way 
mutual compensations [88,41,11]. An up-to-date survey of the current research in 
these directions can be found in ref. [42]. 

The difficulties listed above make the series over periodic orbit conditionally con- 
vergent at best. Nevertheless the existence of a number of special examples where 
the trace formula has been successfully applied encourages to take very seriously 
the insight it offers in the quantum behaviour of classically chaotic systems. In the 
words of Gutzwiller, "as physicists, we have to make a compromise between logic 
and intuition" [87]. 
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2 Quadratic path integrals and functional determinants 

Most systems of physical relevance are described by Lagrangians with strictly pos- 
itive definite kinetic energy. Under such an assumption the quantum propagator 
can be written as a configuration space path integral. The semiclassical approx- 
imation of quantum observables requires the explicit evaluation of quadratic path 
integrals. In configuration space the task is equivalent to computing the determinant 
and the index of an elliptic second order linear differential operator. The index is 
here defined as the number of negative eigenvalues. Forman's theorem reduces the 
computation of the functional determinant to that of the determinant of the funda- 
mental solution (Poisson map) of the linear homogeneous problem associated with 
the nuUspace of the elliptic operator. The theorem is based on the construction of 
homotopy transformations between elliptic operators. Using the same positive def- 
initeness assumption, Morse theory permits to classify different homotopy classes 
and in this way to compute the index. 



2. 1 Path integrals and Lagrangians 

Feynman [62,64] introduced path integrals as fundamental objects governing quan- 
tum dynamics. In non relativistic quantum mechanics Feynman's intuition has now 
evolved into a rigorous mathematical theory [1,47,48,49,50,51,121,32]. 
In the present work quantum mechanical path integrals are derived from an ana- 
lytic continuation of the Wiener measure. The use of the analytic continuation has 
the advantage to unravel the probabilistic interpretation of Feynman path integrals. 
Furthermore it provides a unified formalism for quantum and statistical mechanics. 

The connection between the Feynman path integral and the Wiener measure can 
be established starting form partial differential equations. Consider a continuous, 
forward-in-time stochastic process with support on the same configuration space 
9Jt where the Schrodinger equation (3) was defined. The fundamental object de- 
scribing the stochastic process is the conditional probability also called transition 
probability density Kz to find the process in a point at time T given the position 
at a previous time T'. By definition K-^ transforms as a scalar with respect to the 
invariant measure of SOT under a change of variables Q ^ Q 



I d'Q^I^)K,{Q, T I g', T') = f d''Q^fj{^)k,{Q, T \ Q', T') (22) 

In the presence of a drift field and of a damping potential 0, the transition prob- 
ability is governed by the covariant Fokker-Planck equation 
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dK, z 



dT h 



K,, V T > T' > oo 



2 c M ' 

_ m z 

lim KM T I Q\ T') = ^^'^f'^J^'^ (23) 

The link between the Fokker-Planck equation and the Schrodinger equation is es- 
tablished by an analytical continuation in the parameter z usually referred to as 

Wick rotation. If z is rotated along the unit circle of the complex plane from the 
real to the imaginary positive semi-axis, the identifications 



m 

= C/+^ + ^V„^« (24) 

recover the Hamilton operator (3). In the second of the (24) the momentum operator 
has been replaced with a covariant derivative in order to emphasise that it acts only 
on the vector potential A"'. 

A physically convenient picture of the Fokker-Planck equation is achieved by ap- 
plying to it with the method of the characteristics. The second order spatial deriva- 
tives in (23) impose the characteristics to be solutions of stochastic differential 
equations [49,96,124,91]. If the metric is time independent, it is possible to write 
covariant equations for the characteristics of (23) in the guise of the system of 
Stratonovich stochastic differential equations: 



rfg"(t) = v'-{q{t),t)dt+^^at{t)odwu{t), q-{T') = Q'- 

d^nt) = -T-Mt))<it)odq^{t), g-P{Q') = {a^4){r) 

(25) 

where q°' are the coordinates of the position process while q describes the damping. 
The {(T^lfc^i form a set of vielbeins parallel transported along a path q°'{t) by the 
Christoffel symbols F^^ specified by the metric Qap (appendix A). The vielbeins 
project on OJl the increments dw'^ of a Wiener process (Brownian motion) based on 



(wfc(t)) = Vt < 

(dwkit) dwi{t')) = 5ki5{t- t') dt (26) 
Finally the symbol o highlights the Stratonovich's mid-point rule: 
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aUt) o dwk := lim ^kit + dt) + a^t) ^^^^^ ^ _ ^^^^^^ ^^7) 

dt 10 Z 

The relevant feature of the mid-point discretisation is to cancel terms of the order 
0{dwl) ~ 0{dt) in time differentials [49,91,96,124,149]. Ordinary differential 
calculus therefore applies to Stratonovich stochastic differential equations. More 
details on the geometric meaning of (25) can be found in [96] and appendix B. Here 
it is enough to stress that on a Riemann manifold Wl covariant characteristics for 
the Fokker-Planck equation (22) can be written only by means of path-dependent 
vielbeins. 

Thinking in terms of characteristic curves evinces the intuitive content of the Feynman- 
Kac formula [49,91,96,124]. The transition probability density is the average over 
the Wiener measure of the solutions of (25) connecting Q' to Q in the time interval 

[T',r]: 

^l^)K,(Q,T\Q',T') = IVi,{w(t))c;(T)5^^(q(T)-Q) 

= q"{t;T',Q',w{t)) T > t > T' (28) 

The Wiener measure can be thought to have support on continuous paths w in 
[T' , T] with square integrable absolute value. A rigorous and compact discussion 
can be found in [49]. The exact result 



w{T)=W g 2{T-T') 

Vn{w{t)) = J (29) 

motivates the representation of the Wiener measure usually encountered in the 
physical literature [165] 



Vn{w{t)) = V[w{t)]e-JT'2 (30) 

This latter suggests to represent the Feynman-Kac formula directly as a path in- 
tegral over the realisations of the position process. This "change of measure" is 
defined by means of the asymptotic expression of Kz{Q,T\Q',T') for short dis- 
placements of the position process from its initial state. Under reasonable smooth- 
ness assumptions on the drift and the metric the asymptotics is obtained by substi- 
tuting the short time solution of the stochastic system (25) into the Feynman-Kac 
equation and yields 



/ m \ 9 1 rT +dt, „ , 

K,{Q,T'+dt\Q',T')^(-^^) e-^/T' ^^^^M+o{dt) (31) 

\Z7i z natj 
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The Lagrangian appearing in the exponential is [5]: 



C,{q,q) = - q-v^^ + z^(l)+— V,i;" - (32) 
2 2 Dm 



and it is evaluated along the geodesic, supposed unique for dt small enough, con- 
necting Q to Q'. The notation in (32) means 

V«t;" = a„t;- + r-X (33) 

with R the curvature scalar defined by the metric g^fi- The derivation of (31) is 
summarised in appendix B. 

The path integral for finite time differences follows by iterating N times the con- 
volution integral of short-time kernels 

K,{Q, r + 2dt\ Q', T') = {K, ★ K,){Q, r + 2dt\ Q', T') 

:= / d^Q"Jg{Q»)K,{Q, T'+2dt\Q", T'+dt)K,(Q", T'+dt\Q', T') (34) 

and then taking the limit: 

K,{Q,T\Q',T') = lim {K,* ... * K,){Q,T\Q' ,T') 

N ] oo 
N dt=T-T' 

fqiT)=Q . _ , , . ^ fT 

lq{T')=Q' 

The procedure outlined above can be made completely rigorous [5]. 



:/ V[^q{t)]e--^iT''^''^^ (35) 

Jg(T')=0' 



The quantum mechanical propagator corresponds to the analytical continuation of 
(35) obtained by setting z = exp{— i 9} and rotating 9 from zero to 9 equal 7r/2 
[126]. The Feynman path integral for z equal to i is 



K{Q,T\Q\T') = r V[^q{t)]eilT'dtc (35) 

Jq{T')=Q' 

where the Lagrangian 



>C g, q) = ^q-vf -(!)+— Vav'' + -— (37) 
2 2 bm 

with the identifications (24) assumes the form 
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q) = - gapr f + qaA''-U+—- (38) 
2 6m 

In the last two formulae and from now on any reference to z is dropped. The co- 
variant path integral (36), (38) was originally derived by B.S. DeWitt in ref. [46] 
starting from the short time solution of the Schrodinger equation. Note that on a 
curved manifold the covariant Hamilton operator (3) corresponds to a Lagrangian 
comprising a curvature term vanishing in the classical limit (see also discussion in 
[145]). 

For analytical and numerical purposes, it is often convenient to interpret the limit 
of the iteration procedure (35) as the continuum limit of a single N x d lattice 
integral [99] with time mesh 



A direct lattice construction of the quantum propagator is also possible. The La- 
grangian in the exponential is recovered if and 51^/3 are discretised according 

to the mid-point rule [145]. The discretisation rule reflects the interpolation with a 
differentiable curve, a geodesic, which was used to derive the short time transition 
probability density (31). 

A mid-point discretisation permits to import the rules of ordinary calculus for for- 
mal manipulations under path integral sign. The circumstance was early realised 
by Feynman [62] by requiring the gauge invariance of the propagator of a quantum 
particle interacting with an electromagnetic field. 

From the knowledge of the propagator it is possible to reconstruct the dynamics of 
all relevant observables in quantum mechanics. This is generally done by averaging 
over the propagator the kernel of a self -adjoint operator 



(O) = / d''Qd''Q'^Jg{Q)g{Q')0{Q, Q')K{Q, T \ Q\ T') (40) 

If the average operation is reabsorbed in the definition of the functional measure, 
more general path integrals are obtained 



(O) = / V[^q{t)] et Jt' '^^(^+^0) (41) 

To wit, the kernel 0{Q' , Q) will not only modify the potential in the original La- 
grangian (38) but it will also impose new boundary conditions on its lattice dis- 
cretisation. The support of the path integral is then identified with the space of 
continuous paths *P satisfying such boundary conditions. An example of (41) is the 
trace of the propagator: setting 0{Q' , Q) = 5^'^\Q' - Q)/y/g{Q) yields 
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TrK :-- 



(42) 



where L^M is the loop space, the space of quantum paths closed in [T", T]. 



2.2 Semiclassical approximation and quadratic path integrals 



The path integral (42) suggests a physically intuitive picture of the semiclassical 
approximation. The leading order should correspond to quantum paths exploring 
the configuration space Tl only in a neighborhood of classical trajectories satisfy- 
ing the boundary conditions ^. Namely, had the path integral support on the set of 
at least once differentiable vector fields over Wl, 



= {qit) e CW([T', r],mt) I {q{r),q{T)) e ^} 



(43) 



it would be possible to give a literal meaning to the time derivatives in the path inte- 
gral Lagrangian. Once restricted to Cfp, the path integral concentrates for vanishing 
h around those curves qce for which the action is stationary 



= 



6S 



6q-{t) 



Qclit) 



_5g°(t) Jt 



dtC 



The latter condition is equivalent to require qce to satisfy 



(44) 



dC d dC 



dq"' dt 



= 



5?" 



dC 



dq°' 



T 9q^ 



for all fluctuations Sq in the tangent space to qce 

T,,C = {Sqit) e CW([r',T],T,^,2Jl) I {5q{T'),5qiT)) e 5B} 
Thus the representation of quantum paths 



(45) 



(46) 



q''{t)^q^,{t) + VhSq"{t) 
gives the asymptotic "naive" approximation of the path integral on C^p 



(47) 
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(O) ~ e^^T'^oi f V[^Jq(t)]e'lT''^'°le^ + 0{Vh) (48) 

{QceeC} -''^"ct^ 

the sum being extended to all extremal curves in C. According to (48) the path 
integral to evaluate is governed by the second variation Lagrangian 



with Lgg, Lgg, Lgg, d X tlmc dcpcndent real matrices obtained from the second 
derivatives of the original Lagrangian evaluated along the stationary trajectory Qce. 
The arrow over the derivatives indicates that they act to their left. In particular 
the "mass tensor" L^^ is equal to the metric tensor evaluated along the extremal 
trajectory 

(L^^)„^(t) := mg^f^iq^t)) , L^f (t) := iL-l)ap{t) (50) 

and it is therefore symmetric and strictly positive definite. The measure terms ■y/g^ 
are also evaluated along q^t. 

Unfortunately the situation is more complicated. Feynman path integrals as well 
as the Wiener measure concentrate on nowhere differentiable paths. A discussion 
of intuitive appeal of such issue is given in appendix 3 of [34]. A rigorous justi- 
fication of the semiclassical approximation requires more mathematical work (see 
[51,32] and references therein). However, it turns out that formal rules of path inte- 
gral calculus leading to the correct quantum theory can be inferred from the lattice 
representation of path integrals. 

On a finite lattice the stationary phase approximation can be applied to the La- 
grangian (38) or its generalisations discretised according to the mid-point rule. The 
stationary point is seen to correspond to the discrete version of (45). The inte- 
gral over quadratic fluctuations reduces then to a multidimensional Fresnel integral. 
Definition and basic properties of Fresnel integrals are recalled in appendix G. The 
lattice quadratic action is specified hyanNdxNd dimensional symmetric matrix 
Ln. The entries of Ljv are read off the mid point discretisation of the path integral 
action and from the boundary conditions. If the matrix L^v is non singular a straight- 
forward computation performed in appendix G yields 

.W(^) ^ xo, , (51) 
^|DetjvLjv| 

The result requires some explanations. The prefactor is a complex number de- 
pending only on the boundary conditions 05 imposed on the lattice. 



d d d I t ^ 



qq 



6q^ 



(49) 
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The symbol Dctjv indicates a redefinition of the determinant. Namely, on a finite 
lattice with mesh (39) the propagator path integral gets an overall normalisation 
constant: 



Hi 



N 



-i IL-\ 

me 2 



2 



_2ndt h 

The DetAT operation is defined by setting 



(52) 



^ I ^^^^^ ^^-3^ 

DetivLiv oc „ (53) 

the proportionality factor being N independent and being fixed by >i^. The re- 
definition of the determinant operation is immaterial on a finite lattice but leads to 
great simplifications in the continuum limit. Finally the Fresnel integral is seen to 
produce the phase factor ind^Liv. Due to the normalisation (52) the phase factor 
coincides with the Morse index of L^v: the number of negative eigenvalues of Ljv. 

The lattice computation can formally repeated in the continuum by writing the sec- 
ond variation of the action in terms of a non-singular, self-adjoint Sturm-Liouville 
operator. In general, evaluated on any two smooth vector fields ^, x, the second 
variation defines an infinite dimensional quadratic form 



d d d I t ^ 



qq 



(54) 



a/3 



Integration by parts yields 



<^''5(e,x) = r(Vx)a|?'+ r dtCL^fiX^ (55) 

JT' 



where 



LaPX^ - |(Vx)a + {\-\q)aP ^ + {Uq)apX^ (56) 

and 



(Vx)a {Ui)aP^ +(L4.)a/3X^ (57) 

is the momentum canonically conjugated to x- The notation is motivated in ap- 
pendix E. If the boundary conditions 05 satisfied by the vector fields are such that 
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r(vx)ai?. = o 



(58) 



the differential operation L specifies an operator in L^([T' , T], M'^) self-adjoint 
with respect to the scalar product 

(e,x) jldtC{L^^Ux^ (59) 

The associated eigenvalue problem is 

(60) 

(A„(r),A„(r)) e 05, Vn 

Eigenvalues and eigenvectors in (60) are labelled by quantum numbers collectively 
denoted by n. The operator L<b inherits from the kinetic energy of the Lagrangian a 
strictly positive prefactor of the highest derivative. The fact together with physically 
reasonable smoothness assumptions, permits to apply known results in functional 
analysis [52,56] insuring that the eigenvalue problem admits a countable number of 
solutions orthonormal with respect to (59) with eigenvalues bounded from below 
and accumulating to infinity. The Morse index ind~L<B is therefore well defined. 
The condition (58) is satisfied by the boundary conditions encountered in classi- 
cal variational problems encompassed by Morse theory. Examples are Dirichlet, 
periodic, anti-periodic or focal boundary conditions. The unitary evolution laws 
of quantum mechanics renders these examples the ones typically encountered in 
semiclassical asymptotics. 

The continuum limit of a quadratic path integral can be thought to have support on 
the space of square integrable fluctuations: 

L^([T',T],M^) : = 

{dqit) e [T', T]xR''\ {dq, dq) < oo , { dq{T'), dq{T) ) e } (61) 

Any element of [T" , T] , R'^) can be represented as a series over the eigenvectors 



Sq^'it) = Y.{>^n,Sq)X^{t) (62) 

n 

and therefore 



5''S{5q, 5q) = {5q , L^5q) = {L^5q , 5q) (63) 

The path integral measure consists then of the countable product of Fresnel inte- 
grals over the amplitudes {Xn ,5q). The normalisation of each integral can be in- 
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ferred from the lattice regularisation, compare with appendix G and H. 
The above picture of path integration requires L<s to be non singular: 



KerL<B = (64) 

The hypothesis is too restrictive for most physical applications and in chapter 4 it 
will be shown how it can be loosened. For the rest of the present chapter, however, it 
is supposed to hold. Self-adjoint Sturm-Liouville operators governing the evolution 
of quantum paths in the neighborhood of a classical trajectory are often referred to 
as fluctuation operators. 

In conclusion the leading order of the semiclassical approximation reads 



{6)^ ^ e^T'^c. / V[^Jq{t)]e'lT'^'^le^ (65) 

The result is exact in the case of quadratic theories. It remains therefore the problem 
to evaluate explicitly 



g-«f ind 



iL|([T',T],Rd) ^ ./|DetL<B| 



This can be done by combining classical Morse theory with Forman's theory of 
functional determinants which will be expounded below. 



23 Functional determinants: from physics to mathematics 



The functional determinant of an infinite dimensional operator cannot be defined as 

the product of the eigenvalues. These latter are unbounded and the product would in 
general diverge. Fortunately the physical motivation to introduce functional deter- 
minants helps to overcome the difficulty. To streamline the notation, tensor indices 
will be from now on omitted whenever no ambiguity can arise. 

Functional determinants come about as a result of the infinite dimensional Gaussian 
or Fresnel integrals entailed by quadratic path integrals . Equation (51) above, shows 
that Det n is the quantity the continuum limit of which is really needed. Det n has 
chance to converge in the continuum limit as the result of the mutual cancellation 
of two infinities. This is indeed the case each time it is possible to make sense of 
path integrals according to the definition given in (35). It is therefore appropriate to 
identify the functional determinant Det with the continuum limit of Det at. 
It is a remarkable fact that in all cases when functional determinants exist according 
to the definition just given their expression can be as well recovered if another. 
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mathematically more satisfactory, definition is adopted. Such definition is supplied 
by the ^-function determinant of Ray and Singer [134]. 

Consider a non singular linear (partial) differential operator O which is elliptic on 
a closed domain of definition. The attribute elliptic basically means here that the 
scalar or tensor coefficient of the highest derivative is strictly different from zero. 
The zeta function of O is the series 



iJ":^TrO-', Res > (67) 

ij e SpO 

summing the eigenvalues of O. The symbol Tr in bold characters betokens the ab- 
stract trace operation on O. The series is convergent for Re s large enough. Since a 
finite dimensional square matrix with non-zero eigenvalues satisfies the equality 



IndetO = - ^TrO"^ 
as 



(68) 



s=0 



with now Tr the finite dimensional trace operation, it is tempting to surmise 



(69) 

s=0 

Although Tr does not converge for Re s near zero, the definition is consistent. 
Theorems of Seely [146] insure that, under general conditions, Tr can be ana- 
lytically continued to a meromorphic function of the entire complex s-plane which 
is holomorphic in the origin. 

The analytic continuation is most conveniently achieved by means of the Mellin 
transform of the fundamental solution of 



InDetO := - ^TrO" 
as 



{du + O)f = 

f\u=o = I (70) 
If O is reasonably smooth 



oo 



0-':=-—l duu'-'e-^'' (71) 
is well defined and therefore 

1 f-"^ 



TrO" 



/■oo 

/ duu'-^ Tre-''^ ill) 
Jo 



r{s) 

can be used to define functional determinants 
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The definition through the ^-function has the advantage to provide with a precise 
mathematical meaning the manipulations necessary to compute functional deter- 
minants without diagonalising the operator. The idea is to write functional deter- 
minants in terms of quantities of lower functional dimensionality. This is possible 
if the elliptic operator O is thought as a member of a one parameter, here denoted 
by r, family of elliptic operators Or acting on the same functional domain. As r 
varies from, say, tq to ri the elements of the family are smoothly deformed from 
to O-,-^. Along the deformation, functional determinants satisfy the differential 
equation 



drhiDet Or = Tr[{drOr)Or^} (73) 

introduced by Feynman in [63]. On the right hand side the differential operator 
drOr acts on the Green function O^^. The differential equation translates into a 
useful prescription to compute Dct O if the initial condition Dct Oq is known. This 
is never a problem because there are examples of exactly solvable quadratic path 
integrals from which an initial condition can be read off. There are instead two 
other difficulties with the functional determinant flow equation (73). 
The first problem is related to the exact meaning of the trace operation. The discus- 
sion can be made more concrete for the family of Sturm-Liouville operators acting 
as 



Lr)'^ '■ — 

_ Z I "Tf' i _|_ I _|t ^ ^_ I I "T ('I _ 

dt'^ 99;-^ ^qq;TJlP '-qq \'-qq;T 



L5,;r)7/3(74) 



with 



Vt e [T',T] 



(75) 



on square integrable vector fields satisfying some boundary conditions OS in [T', T]. 
Note that the tensor prefactor of the second order derivative in (74) is r indepen- 
dent. At variance with the previous section, it is convenient to loosen for a while 
the hypothesis of self-adjointness. A broader class of Sturm-Liouville operators 
is encompassed if the boundary conditions are described by a pair, not necessarily 
unique of 2 d x 2d matrices (3^i , 3^2) such that the Green functions satisfy 



yi 



+ 3^2 



The simplest example is represented by Cauchy boundary conditions €: 



(76) 
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= (77) 

Cauchy boundary conditions will prove useful in what follows. 

By inspection of (76) the derivative of the Green function displays a discontinuity 

on the diagonal at t equal t'. In Cauchy 's case in particular 

!,,„^l;'(M') = -i. (78) 

An elliptic operator is said to be trace class when the associated singularity does not 
affect the determinant flow. Concretely Sturm-Liouville operators //<8,r are trace 
class when Lqq, \-qq-T are equal to zero. Physically this is the case of quadratic 
fluctuations of a kinetic plus potential Lagrangian in flat space. However, if the 
metric becomes nontrivial as in non Cartesian coordinates or if a vector potential 
sets in, Feynman's equation in the above form is ill-defined. In the physical litera- 
ture the problem is usually solved by coming back to the the lattice regularisation 
of the path integral and showing that the singular terms dash out in the continuum 
limit, see for example [13]. The cancellation can be seen directly on the continuum 
in a more general, boundary condition independent framework due to R. Forman 
[67]. He proved that for a family of elliptic ordinary differential operators of n-th 
order the correct definition of the trace from the analytical continuation (71) is 

dr In DetO^ = - hm Tr drO'' = 

IjtTr {2{t) l\mdrOAt)0^],{t,t') + [l-2{t)] \\mdrOr(t)0^],(t,t')^ 
\\d if n is even 

Xt) = \l (79) 
if n is odd 

The projector 1 weights the contribution to the trace of the limits from above and 
below. For configuration space fluctuation operators L^., n is equal to two. For gen- 
eral odd order differential operators, the projector 2 has a complicated expression. 
Examples of physical relevance of odd order differential operators are Dirac ones 
with n equal one describing quadratic fluctuations in phase space [135,165]. For- 
tunately, simplifications occur in that case since 2 reduces to zero or the identity 
when the matrix prefactor of the highest derivative in O is the identity. The proof 
of (79) is rather technical, the interested reader is referred to [67] for details. 
Cauchy boundary conditions provide the simplest application of Feynman's deter- 
minant flow equation (74). The family of Sturm-Liouville operators (74) is seen to 
share the same functional determinant if Cauchy boundary conditions are imposed 
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The result will be exploited to construct explicitly functional determinants with 
general boundary conditions. 

The second difficulty arising from (73) concerns the classification of the admissible 
smooth deformations or, in the mathematical jargon, homotopy transformations. 
Although the problem is general the discussion will focus on the Sturm-Liouville 
family (74) with self-adjoint boundary conditions 23. The smooth dependence of 
the spectrum of L(b,t on r may produce a divergence of the Green function if a 
zero eigenvalue is crossed during the deformation. Divergences are avoided if the 
homotopy connects operators having the same Morse index ind~ L<s. The loos- 
ening of the condition leads to classify the homotopy transformations by means 
of the Morse indices of the operators they connect. Namely, the Morse index can 
be interpreted as the winding number of the phase of the functional determinant 
DetLsB -r. Two functional determinants are then said to pertain to the same homo- 
topy class if they have the same winding number. On the other hand, if in principle 
the knowledge of the index is required in order to construct functional determinants 
from (73) it is also possible to proceed in the opposite direction. In chapter 3 it will 
be shown that if a self-adjoint Sturm-Liouville operator L>b is smoothly deformed 
by the introduction of a positive definite, self-adjoint perturbation, its Morse index 
can only decrease or remain constant. Hence the number of negative eigenvalues is 
also specified by number of zeroes, counted with their degeneration, encountered 
by the determinant as the coupling constant of the positive definite perturbation is 
increased until the overall Sturm-Liouville operator becomes positive definite. 
For the moment the knowledge of the index will be assumed and the attention 
focused on the derivation of the explicit expression of functional determinants of 
Sturm-Liouville operators. 



2.4 Forman's identity 



Forman's identity [67] relates the ratio of the functional determinants of an elliptic 
operator O with boundary conditions 05 1 and 052 to the determinant of the Pois- 
son map of O. The Poisson map is the fundamental solution of the homogeneous 
problem associated to O. Explicitly Forman's identity reads 



^^^^DetO^ = a,lndet^.(^2,^i) (81) 
DetOsBi ,T 

Qr{^2, ^i) is the Poisson map applied to boundary conditions given by the lin- 
ear complement of 05i and projected on the linear complement %2 of the boundary 

conditions 

Although far more general, Forman's identity will be now derived for Sturm-Liouville 
operators. The proof goes along the same lines of the general case with the advan- 
tage that all operations can be given an explicit expression. 
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The functional detemiinant for Cauchy boundary condition permits to write a more 
general determinant flow equation 



dr In 



DctL(j; T- 

T rT 



j I dtdt'5{t-t') Tr {{drLr){t) [L^]^{t,t') - L^^^{t,t')]} 
The above equation is well defined. The difference of two Green functions 



(82) 







(83) 



is solution of the Poisson homogeneous problem L f = and therefore does not 
have singularities on the diagonal. There is a further advantage in writing (82): the 
integral on the right hand side can be performed exactly. The proof requires though 
to reformulate the Green function problem as a first order one in 2 d-dimensions. 
This is possible in consequence of the identity 



J^^ J^^ dtdt' S{t - t') Tr {(a,LO(i) [L-^\At. t') - Ll\{t, t')] } = 
y dtdt' 5{t - t')Tr {{drVr){t) \P^]At. t') - 'Pe^rit: t')] } (84) 



where 



I d 



- ± - ItU L - • + L- - Lt ) 
and V^^^{t, t'), V^\{t, t') are the Green functions in tangent space: 



(85) 



Vrit) V^\{t, t') = 5{t - 1') hd t, t' e [T' , T] 

y,v^\{T',t') + y^v^\{T,t') = 



(86) 



The statement follows by a direct computation. 

Note that according to (86) the boundary conditions OS are "transversal" to the ma- 
trices 3^1, 3^2- These latter describe the linear complement of *B in a sense which 
will be enunciated more precisely in section 3.1. 

The advantage of working in the tangent space is that the equations of the mo- 
tion and the boundary conditions can be treated on the same footing. To wit, the 
homogeneous problem associated to (85) is 
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Vr{t)dAt,t') = 

dr{t' ,t') ^ hd (87) 

Hence, the Green function fulfilling Cauchy boundary condition € (11) is 

v^\{t,t') =-^{t,t')e{t-t') 
1 t > t' 



e{t - 1') 



(88) 

t < t' 



By (87), (88) the Poisson map ^r{t , t') of Vt is seen to obey boundary conditions 
linearly complementary to Cauchy's. 

Green functions satisfying boundary conditions ® are reducible to a linear combi- 
nation of the Cauchy Green function with a particular solution of the homogeneous 
equation 

V^\{t, t') = V^j^t, t') + ^rit, T') m{T\ t') (89) 
since the compatibility condition 



[yi + y2 dr{T , T')] t') ^ -y2 vi],{T, t') m 

admits a unique solution for ^(T', t') each time V<:s^r does not have zero modes. 
Thus the integrand in (84) can be rewritten as 



■n{(s,p,)(«)5,(«,r')*(J',«')} 

= Tv{(a.P,)W5,(,r)^^^3^^^y.P-(r,o) (9i) 

The derivative of Vt along the homotopy can be eliminated using the identity 



dr[Vr{t)^r{t,T')] = (92) 

Furthermore, the circular symmetry of the trace allows to write 

An integration by parts and the ^-function in the integrand of (84) are then used to 
carry over from t to t' the derivatives in Vr- The explicit form of the Cauchy Green 
function (88) substantiates the otherwise formally obvious equality 
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V^]AT. t')Vr{t') = 5{T - t') hd (94) 
Substituting the last equation in (84) one obtains 



Therefore, provided the homotopy transformation does not encounter zero modes, 
the Forman's identity (81) becomes: 



2.5 Applications of the Forman identity 



The physical content of Forman's identity is to assert that the leading contribution 
of quantum fluctuations to the semiclassical expansion is determined by the lin- 
earised classical dynamics around the stationary trajectory. Canonical covariance 
of classical mechanics is better displayed in phase space. Thus, it is convenient to 
rephrase Forman's identity in terms of phase space quantities. 

The linearised tangent space flow ^(T , T') is connected by a linear transformation 
to the linearised phase space flow F(T, T'): 



(97) 



where the matrix %{t) is defined in appendix E.1.3. The linear phase space flow 
F(T, T') is specified by the Poisson brackets (appendix E) of the flow solution of 
the first variation 



F(r,r') 



{gc£(r),Pe£(T')}p.b. 

{Pct{T) ,l?c£(T')}p.b. 



-{Pc£(T) ,gc£(T')}p.b. 



(98) 



where and pd are classical positions and momenta and each bracket symboli- 
cally represents ad x c? set of relations. 



In terms of phase space quantities the Sturm-Liouville problem is restated as 
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T' <t < T 



Yi 


A„(r') 


+ Y2 


A„(r) 


= 


(99) 




(VA„)(T') 




(VA„)(r) 







while Forman's identity becomes 

DetL<B([r', T]) = x<Bdet[Yi + Y2 F(r , T')] 



(100) 



If the homotopy transformation takes place in the space of self-adjoint operators, 
the prefactor xsg can be evaluated on any reference operator Lsg 



det[Yi+Y2F,„(T,T')] 



(101) 



by keeping track of the change of the Morse index along the homotopy path. The 
situation is illustrated by some examples. 



2. 5. 1 Dirichlet boundary conditions 

Dirichlet boundary conditions are the most common in applications. They are obeyed 
by fluctuations along a classical trajectory connecting in \T' , T] two assigned po- 
sitions in space. Therefore they are associated to the semiclassical approximation 
of the propagator. In the Dirichlet case the phase space matrix pair in (99) can be 
chosen as 



Yi 





















(102) 



The functional determinant is therefore 



DetLDi..([2^', T]) = XDi,. det[- {q,,{T) , ge£(T')}p.b.] (103) 

If the Lagrangian £ describes a single particle with mass m in flat space the pro- 
portionality constant is most conveniently evaluated from free particle motion. The 
formulae given in appendix H give in one dimension 



Det 



-m 



Dir. 



-{?S(r),gf,(T')}p.b 



{[T',T])^ 



27Th{T-T') 



m 



e'2 , 



m 



(104) 
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Quantum fluctuations for free particle motion are quadratic with positive definite 
spectrum. Therefore the Morse index is equal to zero. When such information is 
combined with (104) one finally gets into 



DetLDi..([T',T]) = 

(2 TT ny I det{g,,(r) , gc^(r')}p.b.| (105) 

The result is usually presented in a slightly different form. The action integral eval- 
uated along becomes a function of the initial and final point of the trajectory 



S{Q,T,Q',r)= [ dtCiq,e,qce,t) (106) 
Jt' 

The action with the above functional dependence is often referred to as Hamilton 
principal function. The partial derivatives of S{Q, T; Q', T') are proven to fulfill 
[7,98,123] 



dS 



a 



dS 

PL Pee AT') = (107) 

The total differential of the initial momenta with respect to the positions at time 
interval ends is obtained from the second derivatives of the Hamilton principal 
function 



dp. AT) = - "Q'g^i^. - "V'scMr- 

The Poisson brackets {QciiT) , Qcei'P')}p.h. yield the displacement of the final posi- 
tion QceiT) versus a change of the initial momentum pc^/3(7") while all other initial 
momenta and positions are kept fixed. This is exactly the same as the inverse of the 
differential (108) when the Q'^'s are held constant. It follows 



det- 



^ (109) 



|det[{ge£(T),ge^(T')}p.b 



Thus the familiar form of the semiclassical propagator is recovered 



K{Q,T\Q,T') 



det 



dQf^dQ"^ 



,i K [ind-Loir.{[T',T])+i] 



[2n7r]'^/2 



(110) 
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The semiclassical propagator is named after Van Vleck [155], Pauli [127] and 
DeWitt-Morette [117] who rediscovered it in independent contributions. The re- 
lation with functional determinants with Dirichlet boundary conditions was first 
proven by Gel'fand and Yaglom [72]. Finally the complete form of the semiclassi- 
cal approximation featuring the Morse index for open extremals is due to Gutzwiller. 

The expression of the semiclassical propagator in curved spaces is also known 
[46,102]. In order to recover it from Forman's theorem one needs to evaluate the 
functional determinant in a reference case. The free particle supplies again the 
needed information although more work is demanded in comparison to the flat 
space counterpart. The interested reader can find the details of the computation 
in appendix 9.2 of ref. [99]. It turns out that the resulting functional determinant 
predicts the semiclassical asymptotics 



K{Q,T\qf,T') 



det 



• S(Q,T,Q',t') 



■if [ind-LDir.([T^'>T])+f] 



\g{Q)\H^f^^Y'^\9{Q')\^ 



(111) 



The result coincides with De Witt's computation of the quantum propagator [46] 
holding for arbitrary time separations. For small times the Morse index is zero as 
will be discussed in chapter 3. Moreover for small times and small spatial separa- 
tions [50,145] the equality holds 



T -T 



d 



m 



1 + 



det 





1 


dQl^dQ'^ 
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{Q - Q'YiQ - Q'Y + o{AQ^) 



(112) 



in the limit 



0{dT) ~ 0{Q - Q'Y i 



(113) 



Thus the small time expansion of section 2. 1 is recovered if S is identified with the 
action specified by the classical part, h set to zero, of the Lagrangian (38). 



2.5.2 Periodic Boundary conditions 

Periodic boundary conditions are those of interest for trace formulae. They are 
enforced by the choice 



Yi = -Y2 = I 



2d 



(114) 
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Forman's identity yields 



DetLper.([T', T]) = xper. det [ha - F(T, T')] (115) 

If the Sturm-Liouville operator arises from the second variation evaluated on a 
prime periodic trajectory, the time difference = T — T' is the period of the 
orbit. The attribute prime means that the trajectory has wended its path along the 
orbit only once. The matrix 



M := F(T,T') (116) 

is referred to as the monodromy matrix. The stability of any trajectory on the orbit 
is governed by the same M [95,71,1 16,7,123]. A short summary of the properties of 
the monodromy is provided in appendix E.2.2. As it will be shown in chapter 4, the 
monodromy matrix acquires an eigenvalue and a generalised eigenvalue equal to 
one in correspondence of any one parameter continuous symmetry of the classical 
trajectory. The functional determinant (115) is therefore zero. Forman's identity 
can be applied generically only to Sturm-Liouville operators stemming from non- 
autonomous Lagrangian. 

The argument does not apply to oscillators. In a generic time interval [T', T], stable 
oscillators do not admit classical periodic orbits. In consequence the functional 
determinant is non zero for it receives contribution only by quantum fluctuations. 
The Lagrangian of a one dimensional unstable oscillator is also a positive definite 
quadratic form 



^ ^ ^ muj^ ^^^^^ 

so the Morse index is equal to zero. The unstable oscillator provides the reference 
case in flat spaces wherefrom xper. can be determined. By comparison with the 
path integral formulae of appendix H: 



Det 



d^+u;\^J[T',T])^Asmh^u 



rp_rpr 



det [I2 - F(r, T')] = - 4 sinh^ (u 




it is found 
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2.5.3 Smooth deformation of the boundary conditions 

A final example of applications of the Forman identity is the variation of functional 
determinants under rotations of the boundary conditions [67]. 

Under a rotation of phase space such that 

x^^\t) ^ R{t) x{t) (119) 
the linear flow transforms as 

F(^)(T,T') = R-\T) F(T,T') R(T') 
Accordingly, the functional determinant of (99) obeys 

DetL^\[r,T]) _ det 
DetL^{[T',T]) ~ ~det [Yi + Y2 F(T, T')] 



Yi+Y2F('^)(r,r') 
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3 Intersection forms and Morse index theorems 

The construction of functional determinants by means of Feynman's equation (73) 
requires an homotopy classification of the determinant bundle. In the case of self- 
adjoint Sturm-Liouville operators the Morse index provides the classification. 

Morse indices are topological invariants. The circumstance gives great freedom in 
their evaluation. The classical index theory devised by Morse [118,112] provides 
the explicit expressions of the indices associated to general self-adjoint boundary 
conditions. The "Sturm intersection theory" introduced by Bott [29] and further 
refined by Duistermaat [56] and Salamon and Robbin [136] evinces the equivalence 
of Morse theory with the topological characterisation (Maslov index theory) of 
the structural stability of solutions of classical variational problems introduced by 
Maslov [109] and Amol'd [8]. 

The "Sturm intersection theory" focuses on the properties of the eigenvalue (Fred- 
holm) flow [137] of self-adjoint Sturm-Liouville operators. The scope of this chap- 
ter is to provide an elementary review of the intersection theory. The point of view 
is complementary to one adopted in the literature on the Gutzwiller trace formula 
to (see for example [105,106,39,140]) where Maslov indices are deduced from the 
properties of metaplectic operators. 

3.1 Morse index and symplectic geometry 

Let L|( [T' , T] , E'') be the space of square integrable vector fields in x [T , T] 
satisfying some assigned boundary conditions 03. The differential operation 

^ ~^d-^ - \-qq{ \-qq + \-qq;T - \-qq-T)-^ + ^qq {\-qq;T - ^-qqir) (121) 

acting on L^([T' , T], M*^) defines a Sturm-Liouville operator L<b. The operator 
is self-adjoint with respect to the scalar product (59) if for all ^, x belonging to 
L^([T' , T], M*^) the boundary form on the right hand side of the equality 

(X, = {Lix , - (X^ VO I?, + i^Vx) I?. (122) 

vanishes. All Sturm-Liouville operators Ltg considered in this chapter will satisfy 
the condition. 

In ref. [29] R. Bott restated the requirement of vanishing boundary form in a way 
explicitly invariant under linear canonical transformations. He observed that the 
boundary form can be recast as the difference of skew products of the phase space 
lifts of the vector fields ^, x- 
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-(x^VOW + (e^Vx)W = (A^oxJAoOW, t^T',T (123) 

where J is the symplectic pseudo-metric in Darboux coordinates (appendix E), 



-\d 
\d 



(124) 



and A denotes the lift operation. This latter associates to configuration space vector 
fields their canonical momenta. For a linear theory it means 



Aox 



X 

Vx 



(125) 



Hence L<b is self-adjoint if for all ^, x in [T' , T] , W^) 

= (At o xJ A o e)(r) - (At o xJ A o O(r') (126) 
The requirement is fulfilled if the boundary conditions *B admit the representation 



Aox(r') = Ziy 

Aox(T) = Z2y 



Zt J Zi = Zt J Z2 , Zi , Z2 e R^'^ X R2d(i27) 



for some y in W''^ and matrices (Zi , Z2) independent of the interval end times T' , 
T. In the literature the matrices (Zi , Z2) are often referred to as Bott's Hermitian 
pair. 

Bott pairs characterise geometrically self-adjoint boundary conditions ^. 
A Sturm-Liouville problem in d dimensions is completely specified by 2d linearly 
independent boundary conditions. If they are self-adjoint, they define in R^*^ the 
2d-dimensional subspace 



D4rf 



Zi X 
Z2X 



, Zt JZi = Zt JZ2 , Zi, Z2 e M^'^ X R^'^ } (128) 



Furthermore, in R'^"' the condition (126) is equivalent to the vanishing of the non 
degenerate symplectic form "2 obtained by evaluating J © (—J) on any two elements 
of 05 
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A o x(T') 

Aox(r) 



J 
-J 



A ° e(T') 
Aoe(T) 



(129) 



The maximal number of linear independent vectors annihilating a non degenerate 

symplectic form like (129) is exactly 2d in a 4(i-dimensional space [7]. A manifold 
is called Lagrangian if it contains the maximal number of linear independent vec- 
tors pairwise annihilating a symplectic form. Hence 03 is a Lagrangian manifold, 
specifically an hyper-plane, of W^^ with respect to '2. 

In the enunciation of Forman's theorem, the boundary conditions associated to Lsg 
were imposed in the form 



YiAoA„(r')+Y2AoA„(r)=0 



(130) 



for any eigenvector A„ of L^q. The geometrical meaning is a transversality condition 

in \ 





t 




Ylx_ 




Z2I/ 



0, 



t,2d 



(131) 



In other words, the matrix pair (YJ , Y|) describes a linear complement of the bound- 
ary conditions 

Also the zero modes of Lsb admit a symplectic geometry characterisation. The sym- 
plectic form "2 vanishes on vectors framed by the graph 



GrSy := 



y 

Sy 



h2d 



(132) 



of any symplectic matrix S 



S^JS = J (133) 

Any flow F(t, T') solving a linear Hamiltonian system draws a curve in the sym- 
plectic group Sp{2d) (appendix E). In particular f{t,T') can be identified as the 
phase space flow associated to the 2d linear independent solutions j{t) of the ho- 
mogeneous problem 



{Ljnt) = 

fit) = F^(t,T')/(T') + f'i{t,T')m^{T') , a, (3^1, ...,d (134) 
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In variational problems the j(i)'s are called Jacobi fields. The existence of zero 
modes of Ltg corresponds to a non empty intersection between the Lagrangian 
manifolds and GrF(T, T'). By Forman's theorem and the transversality relation 
(131), intersections are analytically characterised by the equivalence 

03 n GrF(T, T') 7^ ^ det[Yi + Ya F(T, T')] = (135) 

The description of zero modes as intersection of Lagrangian manifolds plays a 
major role in the explicit construction of Morse indices. 

The semiclassical approximation brings about Sturm-Liouville operators such that 

5'5(x,0 = {x.L^i) (136) 

holds true. The equality offers a direct way to compute the Morse index of 
[29,56]. The bilinear form can be rewritten using canonical momenta 

<^'«5(x,0 = j^dt [(Vx)^L7ive + X^(L,,-L5,LTiL,,)e] (137) 
The replacement 

U«^L,9 + -^U (138) 
renders the quadratic form 

5'S,{i,i) = £dt [(V0^LTive + e^(L,, + rU-Lt^L7iL^je] (139) 

positive definite for values of the coupling constant r large enough provided U 
and Lqq are strictly positive definite. If the boundary conditions 03 are non-local 
U should be chosen such that L^b ,t remains well defined. For example, if Lsb is 
T — T"-periodic then any matrix time dependent matrix U such that 

U(t + r-r') = U(t) (140) 
is admissible. 

The Morse index receives contribution each time the deformation path encoun- 
ters a Sturm-Liouville L^^^^ with zero modes. Since an arbitrary vector field ^ in 
[T' , T] , W^) is independent of r 

drS'Sr{^,0 = (e,UO > (141) 
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and the index decreases monotonically for increasing r. 

It is instructive to see how (141) comes about by looking directly at the eigenvectors 
of L<s. For the sake of the notation, the parametric dependence on r at r equal to 
zero is from now on omitted. The derivative of an eigenvalue of L^q^t- reads 



dAn,r) (142) 

where An ,t denotes the eigenvector 



(i^<B,rA„,.)(i)=4,rAn,rW, (A o A„ ..(T') , A o A„ ..(T)) gQS (143) 
A double integration by parts in the last addend in (142) gives 



dr{K,r,K,r) + {Xn ,t , ^ K ,r) + "^{dr K ,r , K ,r) (144) 

The first term is zero because eigenvectors are normalised to one at every r, the 
third term vanishes because for all r the eigenvectors satisfy the boundary condi- 
tions. The final result is familiar in quantum mechanics: 



drL,r = {K,r,^K,r) (145) 

Thus the Morse index of L^b coincides with the number of zero modes encountered 
during the deformation for r larger than zero; the contribution at r equal to zero 
would add to the Morse index the nullity of LtQ. 

In practice the prescription requires to solve for arbitrary r the linear Hamiltonian 
system 



Fr{T',T')^\2d (146) 

with 



Kit) = 



-qq 



+Lt IT^I 

^ iq QQ QQ 

^qq ^QQ 



rU 



QQ QQ 



it) 



(147) 



The flow generates the family of linear Lagrangian manifolds Gr Ft-(T, T') of R*^*^. 
The intersections with *B ~ (Zi , Z2) yield the index: 
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ind-L<B([T', T]) = Yl nul {Z2 - F,(r, r')Zi} 

{T>0|'Bn Gr F-r {T,T')} 

^ nul{Yi + Y2F,(r,r')} (148) 

{r>0|'Bn Gr F^(T,r')} 

In the last row the transversality relation (131) has been used. 

The prescription is illustrated by the following elementary example. 



3.1.1 Example: Morse index of the harmonic oscillator 

The Sturm-Liouville operator paired with a one dimensional harmonic oscillator of 
mass m acts on functions in L^([T' , T] , M*^) as 



^2 

= -m — - mu/ (1 - r) 



(149) 



As above, r equal zero corresponds to the Sturm-Liouville operator the Morse in- 
dex whereof is needed. At r equal unity a semi positive definite operator is certainly 
attained. For any r the classical equations of the motion admit the phase space rep- 
resentation 



-1 




(1 




1 




p 





moj^il — T) 





(1 




p 



(150) 



with fundamental solution 



Fr(t,0) 



cos(a;V'l — "Tt) 
muj\/l — T sin(a;-\/l ~ T"t) 



sm(c<j ^Jl—rt) 
vnu3\/\—T 

cos(a;-\/l — 'T'^) 



(151) 



The Morse index can be now computed using Forman's theorem. Consider 

i Dirichlet boundary conditions in [0, T]. 
By (105) the absolute value of the functional determinant is 



|DetLo.,.([0,T])| = 2.ft l^'"'"'^^)l (152) 

The index is given by the number of zeroes for r ranging in ]0, 1]. Zeroes 
occur for 



VT^T^ k7r>0 (153) 
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and are non degenerate. Hence the result of the path integral computation of 
appendix H.0.5 is recovered: 



ind-LDir.([0,T]) =int 



ujT 



TT 



(154) 



The integer part function "int" is defined as the largest integer less or equal to 
the argument, 
ii Periodic boundary conditions in [0, T]. 

The absolute value of the functional determinant is 



|DetLper.([0,T])| = 4 sin^ 



(155) 



and vanishes for 

ljV1-tT = 2A;7r > 



(156) 



The zeroes are doubly degenerate but at r equal to one where the degeneration 
is single. Hence, the result of the computation of appendix H.0.5 is recovered: 



ind"Lper.([0,r]) = 1 + 2int 



u/r 

27. 



(157) 



3.2 Classical Morse index theory 



The relevance of the prescription given in the previous section is mainly conceptual. 
Classical Morse theory [118] permits a more direct evaluation of the index. 
In chapter IV of ref. [118] Morse devises a general formalism encompassing all 
self-adjoint Sturm-Liouville operators. Here instead the attention is restricted to 
self -adjoint boundary conditions such that: 



S''S{x,0 = {x,L^O (158) 

The evaluation of the index is accomplished first for local boundary conditions. 
The paradigm is represented by Dirichlet boundary conditions. The Morse index 
for Dirichlet boundary conditions (open extremals) is derived here by adapting the 
method of ref.'s [103,114]. The result is then exploited to construct the indices 
associated to non local boundary conditions. The procedure is exemplified by the 
periodic case. 



3.2. 1 Dirichlet and local boundary conditions 

Dirichlet boundary conditions are local. Locality enforces self-adjointness in con- 
sequence of the independent vanishing of the two skew products in (122). Therefore 
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any self- adjoint operator Lsq acting on L^([T' , T], W^) can be seen as member of 
a family of operators Ltg ^ parametrised by the length of the time interval [T', r]. 
Consequently, varying r generates a spectral flow. 

For times short enough, Dirichlet boundary conditions render the bilinear form 
(158) positive definite 

{i:L^ir.O\[T',T'+dt] = S'S[T' ,T' +dt]{i: ~ f ) (^') > 0(159) 

as Lgq is strictly positive definite by assumption. 

The spectral flow encounters a zero mode each time a Jacobi field (134) fulfills 
Dirichlet boundary conditions in [T', r]. The direction of the crossing is given by 



T — I JT' 



(160) 



all other derivatives vanishing on a zero mode. The integral localises by means of 
a double integration by parts. The statement is proven by observing that on the 
boundary of the interval [T', r] any eigenvector acquires, beside the parametric, a 
functional dependence on r: 

An,.(T') =0 

A„,.(r) = (161) 
The total differential with respect to r is zero: 

dX 



,T /rril 



dr 



{T') = {drXn,r){T')^0 (162) 



^(^) = + i9rK,r){r) = (163) 

Combined with an integration by parts of (160), equation (163) yields 



< (164) 

since the mass tensor L^^ is strictly positive definite by hypothesis. The Morse 
index is monotonically increasing as a function of the time interval [T', r] length. 

In variational calculus Dirichlet boundary conditions arise from the second varia- 
tion around an open extremal Qce of a classical action S: a trajectory connecting 



r\ n I ( d Xn T I dXn t\ , , 

9.4,.k,.^o = - L,,^J (r) 
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two assigned points in configuration space in the time interval [T, T']. Jacobi fields 
with the property 



j{r) = j(t) = , 



for some r e [T\ T] 



(165) 



are said to have an qce{T')-conjugate point in r. The most famous of the Morse 
index theorems states that: 

Theorem(Morse): Let be an extremum. of S in [T',T] and assume the kinetic 
energy positive definite. The Morse index of 5^»S[r',T] is equal to the number of 
conjugate points to qce{T'). Each conjugate point is counted with its multiplicity. 

The theorem can be rephrased using functional determinants. 

Any linear Hamiltonian flow admits the square block representation 



F{t,T') 



A{tT') B{t,T') 

c{t,r) D{t,r) 



(166) 



(see appendix E.2.2 for details). By Forman's theorem the Morse index with Dirich- 
let boundary conditions is 



ind-LDir.([r',T]) = Yl nulDetLDir.([T',i]) = ^ nu\B(t ,T') (167) 

T' <t<T T' <t<T 



while the nullity is 



nulLDir.([T',T]) = nulB(T,T') 



(168) 



The same construction can be repeated for other local boundary conditions. When 
the spectral flow encounters a zero mode, the Bott representation of the boundary 
conditions (127) permits to write 



dA o X dvn , s 

T' =Zi-f^ T 

or or 



dA o An 
and therefore 



(r) = Z,^(r) 



(a,AoA„,0(T') 
dk o A„ 



dt 



r) + (a,AoAn,0(^) 



(169) 
(170) 



=0 — 



(171) 



en,T=o 
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It might happen that the derivative (171) is zero. In such a case the direction of 
the spectral flow should be inferred using higher order derivatives. The method has 
been thoroughly investigated by Mohring Levit and Smilansky in ref.'s [103,114]. 
However the information carried by the index is of topological nature and does 
not depend on the details of the Sturm-Liouville operator. Thus it is possible to 
rule out degeneration by hypothesis resorting if necessary to small positive definite 
perturbations as in the previous section. In this sense it is not restrictive to consider 
only operators begetting a spectral flow with only regular crossings characterised 
by a non-singular crossing form 



-I(F(i,r'),«):=- Zl (^Ftj^j(i,T')Zi 
= - Z^H(t)Z2 



GrF(t,T')n«B 



GrF(t,T')n<B 



(172) 



In practice, one evaluates the restriction of the matrix on the right hand to the sub- 
space where the crossing occurs. The resulting symmetric matrix has by hypothesis 

non-vanishing eigenvalues. 
The Morse index reads 



ind-L<B([T', T]) = ind- S^S[T',T'+dt] - E sign-I(F(i, T'), 05) (173) 

T'<t<T 
r|GrF(t,T')n<B 

The sign function must be understood as the difference between positive and nega- 
tive eigenvalues whenever multidimensional crossings occur. 
An immediate consequence of (172) is that each time the Bott pair comprises 







(174) 



the Morse index is just the number of crossings counted with their degeneration 



ind-L^([T', T]) = E nul{0 © 1^ - F(t, T')Z^} 



(175) 



T'<t<T 
t|GrF(t,T')e<B 



In such cases the crossing form (172) is determined by the kinetic energy block of 
the Hamiltonian matrix (147) at r equal to zero. 

In practice Dirichlet boundary conditions play a role of preeminence. The reason is 
that Dirichlet fields are (5^<S -orthogonal to any Jacobi field j (134) 



(176) 
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Exploiting such property, Morse indices for general boundary conditions can be 
evaluated proceeding in analogy to section 2.2 where the semiclassical approxima- 
tion was accomplished by separating classical trajectories from quantum fluctua- 
tions. Here the idea is to write any vector field in x in I-'^([T' , T], W^) as a linear 
combination of a Dirichlet ^ and a Jacobi field j 



x{t)=m+j{t) 



(177) 



The Jacobi field is used to match the boundary conditions. If the kinetic energy is 
positive definite the number of negative eigenvalues is finite. Hence, the evaluation 
of the Morse index is reduced to the evaluation of the index of the restriction of 
to a finite dimensional vector space. A general theorem in linear algebra (see 
for example [10] pag. 120) establishes the properties of the index of a symmetric 
form in a finite dimensional real vector space V. In particular for any W C V 
the index is given by 

ind-J^ = ind-J^lw + md-J^\^± + dim{W n W^} - dim{W n KerJ^Xl^S) 

where W"*- := {v e Y\J^{v,w) — ^ w e W}. The theorem formalises the 
obvious observation that the index is the dimension of the subspace of V where the 
quadratic form is 5em/-negative definite minus the intersection of such subspace 
with the kernel of J^. Combined with the decomposition (177) the theorem provides 
the way to evaluate the index of Sturm-Liouville operators with both local and non- 
local boundary conditions. The procedure is illustrated by the periodic case. 

3.2.2 Periodic boundary conditions 
Periodic boundary conditions 

X{t + T-T')^x{t) (179) 

can be considered in [T', T] only if the differential operation L itself is periodic in 
interval [T', T] 



L{t + T- T) = L{t) 



(180) 



A periodic field x admits the representation (177) for ^ a Dirichlet field and j a 
recurrent Jacobi field: 



l2d-F(T,T')) 










.(Vj)(TO_ 




.(vo(r)-(vo(n 



(181) 
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Periodicity can always be enforced in the sense of IJ by a discontinuity of the 
momentum of the Dirichlet fields at one boundary. The sets of recurrent Jacobi 
fields and Dirichlet fields are non intersecting if 



detB(r,T') ^ (182) 

where the d x d-dimensional real matrix B(T, T') is specified by the block repre- 
sentation (166) of the linear flow. Whenever (182) holds true, one has 



ind-Lpe,.(r,r]) =ind- 5'S{j,j) +md-Luir.{[T',T]) (183) 

rec. 

The first term on the right hand is amenable to the more explicit form 

md-S^S{j,j) =ind-(DB-^ + B-^A-B-^-B-^t)(r,r') (184) 



rec. 



The result follows from the recurrence condition (181) and the use of the properties 

of a symplectic matrix given in appendix E.2.1. The same properties insure that 
(184) is the index of a symmetric matrix and therefore it is well defined. Morse 
referred to (184) as the order of concavity of a periodic extremal, see ref. [118] 
pag. 71. 

It remains to analyse the physical meaning of (182). From the point of view of the 
second variation (182) is not a restrictive condition. The stability of a periodic orbit 
coincides with the stability of any of the trajectories covering the orbit. A trajectory 
starting at time T' in generic position will fulfill the condition. 
Non generic situations when (182) does not hold true can treated using the general 
formula (178) with W, W"*- denoting respectively the sets of Dirichlet and recurrent 
Jacobi fields in [T', T], [12]. 

Morse indices appear from the above examples to be determined by the oscillations 
of the Jacobi fields. An intuitive explanation is provided by variational calculus. 
Jacobi fields are extremals of the second variation. Because the kinetic energy is 
positive definite, the second variation as function of the upper time T cannot be- 
come negative on any functional subspace until a Jacobi field does not encounter a 
conjugate point. 

The Morse index for periodic extremals can be also derived by purely symplectic 
geometrical methods based on the analysis of the canonical form of the linearised 
flow [93,94]. 



3.3 Elementary applications 



The simplest applications of Morse index theorems are offered by one dimensional 
Hamiltonian systems. 



41 



2 4 6 8 10 

Fig. 1. In decreasing gray level position, velocity and acceleration of a particle of unit mass 
in the quartic potential V{q) = + The acceleration is minimal (maximal) and 
negative (positive) at odd (even) turning points of the velocity. 

i Particle in a potential well 
Consider the action functional 



S{q,q) ^ [ dt 
Jo 



mq^ 

-Tr-u{q) 



(185) 



The potential U is chosen to be a positive definite analytic function of q. It is also 
supposed to give rise, at least within a certain energy range, to bounded motion. 
The case of a quadratic potentials is ruled out by hypothesis and it will be treated 
apart. A concrete example is provided by the anharmonic oscillator 



f/(9) = ^ + ^ (185) 

Classical trajectories are extremals of the action occurring on the curves of level of 
the Hamiltonian 



E = H{p,q) = ^ + U{q), p= — ^mq (187) 

2 m oq 

which are eventually closed in consequence of the postulated shape of the potential. 
In fig. 1 a typical solution is plotted for the anharmonic quartic potential (186) in 
rescaled adimensional units. 

Let Xce{t; E) be the phase space lift of a trajectory on a periodic orbit of prime 
period Tci{E). A suitable basis for the linearised dynamics is represented by the Ja- 
cobi fields x{t; E), dEx{t; E). Linear independence is guaranteed by the invariance 
of the skew product 
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Fig. 2. The Jacobi fields and dEQci relative to the motion of a particle of unit mass in 
the potential V{q) = + The envelope of the amplitude peaks of the non periodic 
Jacobi field dEQci is linear with slope —dTd/dE. The conservation of the skew product 
imposes the non periodic Jacobi fields to be positive at odd turning points of the velocity. 

dn 



{dExl^Mce){t;E) := {dEPceQce- dEqctPci){t\E) = 



dE 



1 



(188) 



The skew product coincides with the Wronskian determinant of the linearised dy- 
namics. In fig. 2 the two Jacobi fields are plotted for (186). The Jacobi field of the 
second kind dEx{t; E) describes a linear instability [129,51] 



t — T' { dx\ 

dEx{t;E) ■.^K^^x,e{t;E)+ (— j {t;E) 



(189) 



= const. 



with 



dT, 



K :— — 



ci 



dE 



(190) 



The linearised flow in phase space is (compare with formula (E.58) in appendix E) 



F(t, T') = [x,i{t- E), dEx(t; E)] [x,i(T'; E), dEx(T'; E)]'' 



(191) 



The Sturm-Liouville operator associated to the second variation along a trajectory 
of energy E is 



r 



'fir 



{qce{t;E)) 



(192) 



By (191) the Morse index for Dirichlet boundary conditions is determined by the 
number of zeroes of 
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B(t, T') = dEqciit; E) qce(T'; E) - dEqa{T' , E) q,i{t] E) 



(193) 



The right hand side is a linear combination of solutions of the linearised dynamics. 
Hence it is itself a Jacobi field vanishing at t equal T' . The zeroes encountered 

by (193) in ]T, T'[ define the number of conjugate points and therefore the Morse 
index associated to Dirichlet boundary conditions in the interval [T', T]. By (193) 
also follows that if qd starts from a turning point, {qd equal zero), the index will 
be specified by the number of turning points encountered in ]T', T[. 

Let now T be equal to rTd with r a positive integer. Furthermore suppose the 
trajectory q^ to start form a generic point (not a turning point) on the periodic orbit. 
Observe that B(t, T') starts with positive derivative and has the same signature of 
K over multiples of the period. Moreover from the Wronskian condition (188) on a 
turning point 

l^-{dEqctPct){ti.^:,E) (194) 

follows that B(t, T') must change sign at least once over one period, since velocity 
and acceleration change sign twice. Hence 

ind-LDir.(r, T' + rT,,]) = 2r -^-^^^^ (195) 
If instead the trajectory starts from a turning point then 



ind-LDir.([T',T' + rT,,]) =2r-l 



(196) 



Associating periodic boundary conditions to the differential operation (192) in the 
same interval [T', T'+r Tce\ defines the self-adjoint Sturm-Liouville operator Lper.. 
The periodic Morse index is determined from the sum of the Dirichlet index and 
the order of concavity. This latter is read from the monodromy matrix 



M = 



l + '^{E){pdqct){TcuE) 



dE 

^{E)pUTd;E) 



1 - 



dT^i 
dE 

dE 



{E)ql,{T.uE) 
(E) ipdqce)iTcf, E) 



(197) 



If the starting point of the periodic trajectory is taken in generic position, the order 
of concavity is specified by (184) and vanishes identically: 



ind-Lper.([T',T' + rT,,]) = 2r 
nulLper.([T',T' + rr,,]) = l 



1 — signre 



(198) 



If instead the trajectory starts from a turning point, the general index formula (178) 
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Fig. 3. The period T^t of a particle of unit mass in the potential V{q) = /2 + /A 
plotted versus the energy. The period decreases from the value Tc£(0) = 2 7r, the period 
of the harmonic oscillator, with asymptotic behaviour Tci{E) ~ E~ for E tending to 
infinity. 

must be applied. To the index (196) must be added a non zero contribution from 
the order of concavity. The non zero contribution comes from the recurrent Jacobi 
field 



Jrec.(t) = (lci{t- E)dEVci{T'- E) - dEqci^ E)p,,{T'- E) 
with now 



(199) 



ge€(T' + rTe,;£;) = 
and it is equal to 



(200) 



dT 

525(j,ec.,Jrec.) = M21 = -j^{E)pl,{T,,- E) 



(201) 



The sum of the Dirichlet index and of the order of concavity recovers (198). In the 
case of the Hamiltonian (186) the index is seen to be even since the period is a 
monotonically decreasing function of the energy (fig. 3). 

Harmonic oscillator 

The counting of the conjugate points in [0, T] yields (167) 



ind~LDir.([0,r]) = #(zeroes ofsin(a;i) in]0,r[) = int 



TT 



(202) 
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in agreement with the computations of section 3.1.1 and of appendix H.0.5. In 
consequence, the periodic Morse index is 



ind-Lpe,.([0, T]) = ind-LDir.([0, T]) + ind' 



cos{luT) - 1 
sin (a; T) 



1 + 2 int 



2tt 



(203) 



for T generic. 



3.4 Crossing forms and Maslov indices 



The constructions of Morse indices of the previous two sections are reconciled if a 

more general topological framework is developed to investigate the spectral flow. 
This was done by Duistermaat in [56] who refined ideas introduced by Bott in [29]. 

First, the following heuristic argument helps to understand the topological invari- 
ance of the Morse index. The Morse index of a self-adjoint Lsg is the number of 
negative eigenvalues 



ind-L<B= ^i-^j) (204) 

ij e SpL<8 

If L(B is embedded into a smooth n-parameters family of self-adjoint Sturm-Liouville 
operators including a positive definite element, the index can be thought as the line 
integral in the space of the parameters connecting L«b to the positive definite ele- 
ment. Parametrising with r such curve, the index takes the form 



r' 

ind L^=- ^ / drdr9i-£n,r) 

= E r'dT^S{e„,r) (205) 

where it is assumed that zero and r' correspond respectively io and the positive 
definite operator. In order the equality to hold true, only regular crossing of zero 
modes are supposed to occur. The topological invariance of the index stems from 
the identical vanishing of loop integrals over exact differentials which ensures the 
independence of the index of the curve used to compute (205). 

The last equality in (205) states that along any curve in the space of self-adjoint 
Sturm-Liouville operators, the Morse index is specified by the sign of the eigen- 
value flow on zero modes. Proceedings as in section 3.2.1 yields 



46 



9r4,r = (An,r , (^r-^QS ,r) A„ ,r) (206) 

The integral localises on zero modes. Differentiating the Jacobi equation 



{LrJr){t)= (207) 

versus the curve parameter r allows to prove the chain of equalities 

{jr , {drLr)jr) = - {jr , Lrdrjr) = ^(j > ^rjr) (208) 

The restriction to a crossing point defines the crossing form [56,136,137,138]: 

-1(F,(T,T'),S) :=n(j,,a,jOI<BnGrF.(r,T') (209) 

The requirement of smoothness of the deformation can be weakened. Provided the 
spectral flow is continuous, discontinuous non vanishing derivatives of the zero 
eigenvalues also allow to infer the direction of the spectral flow. 
The Morse index of Lsg is 

ind-L<B= E sign-I(F,(T, T'),?B) (210) 

{r>0|53nGrF^(T,T')} 

for any homotopy transformation with only regular crossings connecting Ltg to a 
positive definite Sturm-Liouville operator. The sign of the crossing form over a dis- 
continuous derivative is the arithmetic average of the derivatives from the left and 
the right of the critical value of r. 

The crossing forms (172) produced by deformations of the time interval are par- 
ticular examples. The fact is recognised by looking at the explicit expression of 
crossing forms. A crossing occurs each time for some x,y e R^*^ the equality 
holds 

GrF^x (211) 
(Zi , Z2) being the Bott pair describing *B. Hence (209) is equivalent to 



-i(F.(r,r'),s) = - zl (pt j^j (r,r') 



Zi 



(212) 

GrF^(T,T')n<B 



The matrix on the right hand side is symmetric at glance since smooth deforma- 
tion of symplectic matrices are governed by linear Hamiltonian equations (ap- 
pendix E.2.3). 
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The formulae (211), (212) disclose the possibility to define the Morse index of a 
self-adjoint Sturm-Liouville operator as the index of homotopy transformations of 
finite dimensional quadratic forms. Namely an arbitrary deviation from the crossing 
condition (211) can be written in the form 









+ 











GrF^(r,r')x, x,y,z e 



!>2d 



(213) 



provided the Bott pair (Vi , V2) spans in the sense of section 3.1, some linear com- 
plement QJ in R^'^ of GrF^ and *B. In such a case (213) admits always solution for 



[Z2-F,Zi]y=[V2-F,Vi]z 
On solutions Zr the symplectic form on 55 



(214) 





t 












J 




Ziy 


MiZr 




-J 




Z2y 



defines the family of symmetric matrices in '. 



(215) 



X 



Q(Q3,5J; GrF,) : = 

\ {[ZtjV2 - ZlJVi][V2 - F,Vi]-i[Z2 - F,Zi] + C.C.} 



(216) 



The spectrum of Q(Q5, 2J; Gr F,-) contains a zero if and only if the parameter r at- 
tains a value r* such that Lsg has a zero mode. Furthermore a direct computation 
shows that 



signQ(Q5, QJ; Gr F^*+d^) = signQ(<B, 5J; Gr F^,) + sign~I(F^*, 5B) (217) 

Proceeding in this way, the Morse index of a self-adjoint Sturm-Liouville operator 
is identified with the Maslov index, the index of the graph of a continuous curves 
in Sp{2 d) having only regular crossings with a given Lagrangian subspace OS of 
R.'^'^ [8,90,56,136]. The identification permits to evaluate the Morse index without 
making direct reference to the corresponding operator. In contrast to the eigenvalues 
of the Sturm-Liouville operator, the eigenvalues of the finite dimensional quadratic 
form (216) as functions of r can also change sign by diverging. The circumstance 
insures that arbitrary large values of the Morse index are eventually attained. 

According to Robbin and Salamon [136], the Maslov index of the graph GrS^ of 
any one parameter family of matrices S(r) in Sp{2d) is defined as 
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K(S(r),?B,[ri,r2]) = -signn(S(ri),®) 

+ Yl sign"I(S(T),®) + ^sign-I(GrS(T2),®) (218) 

T1<T<T2 

r|GrS(T)e<B 

End-points contribute of course only if an intersection occurs there. Maslov indices 
(218) enjoy the following properties. 

i Naturality: the index is invariant under a simultaneous symplectic transforma- 
tion * of <B and GrS(T): 

S(r) , * » , [n, T2]) = K(S(t) , 5B , [ti, t^]) (219) 

for all ^ = (^'i , ^^2) leaving invariant the symplectic form H in R^'^ 

*1 J ^1 = J ^2 (220) 

ii Catenation: for ti < T2 < T3 

^(S(t) , ^ , [n, Ts]) = H(S(t) , *B , [n, T2]) + i^(S(r) , ^ , [r„ T3]) (221) 

iii Product: if the curve S(t) = S'^^)(r) © S*^^)(r) and the Lagrangian subspace 

i^(S(i)(^)©S('n^),^^'^ © [n,r2]) = 

H(S(i)(t), S(^), [n, T2]) + ^(S(2)(r), ^(2), [t2, T3]) (222) 

iv Homotopy: the Maslov index is invariant under fixed end-points homotopies. 

Of the above properties the first three are immediate consequence of the definition. 
Homotopy was first proven by Arnol'd in [8]. The proof proceeds by showing that 
the signature of the crossing form defines a Brouwer degree [1 13,68] for the cross- 
ing. The homotopy invariance of the Brouwer degree is then a standard result of 
differential topology [113]. 

One can use Maslov indices to compute the Morse index of a self-adjoint Sturm- 
Liouville operator -L^(a) with respect to the Morse index of a second one L^(6) . 
Under the hypothesis that both *B^"^ and 53 have no intersections with Gr S(t) at 
path end-points ti,T2, Hormander proved in ref. [90] section 3.3 the equality 



^(S(t),®(«) , [ti,T2]) = )k{S{T)M'^ , [Ti,T2]) + 

I fsignQ(Q3(^\ GrS(T2); 05^")) - signQ(Q3(''\ GrS(Ti); S^")) 



(223) 



Replacing S(r) with F(r, T') shows that Hormander's relation is a generalisation 
of the relation (183) between the indices of Dirichlet and periodic boundary condi- 
tions. 
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F(t,T') 

Fig. 4. The homotopy loop 

Homotopy and Hormander's identity (223) provide all the needed tools to under- 
stand the origin of the different prescriptions given in the literature for the evalua- 
tion of Morse indices. The following sections are devoted to illustrate this point. 
The reader interested in a more detailed presentation of Maslov index theory will 
enjoy the very accessible and comprehensive presentation by Robbin and Salamon 
in ref. [136]. The needed basic background in topology can be found in [113]. 



3.5 Topological construction of Morse indices 



Using the results of the previous section, the Morse index of Lsg, 53 ~ (Zi, Z2) can 
be computed from any closed loop formed by varying fr{t, T') in the (r, t)-plane. 
It comprises the following oriented paths 

1 F^(T,T') , dr < T <t' 

The flow ^r{T, T') is solution of (146) with the family of Hamiltonian ma- 
trices generated by adding a positive definite potential term as in section 3.1. 
The initial value dr provides for an arbitrary small positive deformation which 
removes eventual zero modes. The upper value t' is chosen large enough such 
that the Sturm-Liouville operator becomes positive definite. Thus, the Maslov 
index of the path coincides with Morse index of L^: 

«i(F,«B)=ind-L^([r',r])= nul{Z2-F,(r,r')Zi} (224) 

T>0 

<BnGrFr(T,r') 

2 F^/(t,T') , T > t > r 

For T large enough the operator is dominated by the positive definite potential. 
No zero mode can occur until t reaches T' since the quadratic form (139) is 
positive definite for all t strictly larger than T'. At initial time the condition 

Fr{T',T')^hd, Vr (225) 
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may generate an intersection between *B and Gr l2d. The Maslov index (218) 
is in such a case 



K2(F,^)= ^sign{z^H,,(r)Z2 



GrladnOS 



(226) 



since the crossing form can be directly expressed in terms of the Hamilto- 
nian matrix (147). The sign is the opposite than in (209) due to the negative 
orientation of the path. 

3 f{t,T'), r < t < T 

The corresponding Maslov index is 



Grladn® 



- Yl H(t)Z2 } - I sign I Zl H(T)Z2 



T'<t<T 
t|GrF(t,T')n<B 



GrF(r,r')n<B 



(227) 



4 F^(r,r') , < T < 

By construction the last branch may encounter a Lagrangian intersection only 
in the initial point. The Maslov index is therefore 



^^4(F, ®) = - ^ sign I Zt ^Ft J ^ j (T, T') Z, 



T=0 ,GrF(T,r') n <B , 



nul{Z2-F(r,r')Zi} 



(228) 



The last equality holds true because increasing r begets by construction only 
positive definite crossing forms along this branch of the homotopy path. 

The Maslov index of the overall loop is zero 



(229) 



Therefore it follows 



ind-L^([T',r]) = Ki(F,«) = -[K2(F,®) + K3(F,«) + K4(F,«)] (230) 

The Maslov index of the second branch of the loop can be evaluated a priori. For 
any r the Hamiltonian matrix (147) is similar to 



-1 

-qq 



it) 



(231) 
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for a suitable non singular matrix. For r large enough, (231) is the orthogonal 
sum of two d X d dimensional blocks with opposite sign definition. The observa- 
tion recovers the result (175) obtained for local boundary conditions. If Z2 = © 1^ 
the Morse index coincides with the number of crossings counted with their degen- 
eration in the open time interval ]T', T[. Any overlap of the Bott matrix Z2 with 
the position space produces a non trivial generalisation of the order of concavity 
defined in subsection 3.2.2 for the particular case of periodic boundary conditions. 
Homotopy invariance of Maslov indices and Hormander's identity (223) endow 
with great freedom in the choice of the representation of the Morse index which 
most suits the analysis of the physical origin of the contributions. Periodic bound- 
ary conditions well illustrate the general situation. 



3.6 Morse index and structural stability of periodic orbits 



The difference T — T' is in this section supposed to be equal or a multiple of the 
period of the Sturm-Liouville operator. 

The Bott pair describing periodic boundary conditions is (Zi , Z2) = {hd, hd)- 
Therefore the Maslov index of the second branch of the homotopy loop is zero. 
Hormander's identity can be applied directly to the catenation of the second and 
fourth branch. The identification of QS*^"^ with (l2<i, bd) and of ^"^^^ with any 21 ~ 
(Zi, Z2) such that 



2in(i2d, i2d) = 

21 n GrF(r, T') = 



det(Z2 - Zi) 

det[Z2-F(r,r')Zi] ^ 



(232) 



yields 



E^^{F, ihd, l2d)}= ^^3(F,2t) + -signQ(2l,GrF(r,r'); (hd, \2d) ) (233) 

i=2 ^ 

A more detailed proof can be found in theorem 4.3 of ref. [56]. The signature and 
the index of any symmetric form T 'va.2d dimensions are connected by the relation 



ind T — - \ld — Ker.F — sign.F] 
Thus one can recast the periodic index in the form 



(234) 



ind-Lpe,.([r, T'\) = -J^3(F, 2t) + ind-Q(GrF(r, T'), 21; (b,^, \2d) ) - c?(235) 
which is the final expression given by Duistermaat in ref. [56]. 
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The topological expression of the periodic Morse index can be used to analyse 
how the stability properties of the linear flow F(t, T') are reflected in the index. 
The use of homotopy simplifies the analysis. The idea is to embed the linear flow 
F(t, T') into a family of curves in Sp{2d) with fixed end points. By property iv of 
section 3.4 the Maslov index is the same for any element of family. The Floquet 
representation of the linear flow 



F(t,0) = Pe(t) exp|-jy^ dt' H{t')j =: Pe{t) e'-'^ ^ (236) 

is suited for the construction. In order to streamline the notation T' has been set to 
zero. The matrix Pe{t) is supposed to have prime period T, while H is the average 
of the Hamiltonian over a period. The Floquet form is the element of the continuous 
family 



G(t,s) = { LV ; J - - 2 (237) 

Pe(t)Pe-^[s {t - T)] e^s{t-T)+t]^^H ^<t<T 

attained by setting s equal to zero. The end points of the homotopy family are 
independent of s. In order the Floquet representation to be fully specified, it is 
necessary to fix a convention on the phase of the eigenvalues of the monodromy 
matrix: 



M := e^-"^"^ (238) 

A convenient choice is to restrict the phase to the interval [0, 2 7r[. The choice at- 
tributes to the exponential matrix in the Floquet representation a winding number 
equal to zero. Gel'fand and Lidskii introduced in [71] the winding number of a 
linear Hamiltonian flow as a topological characterisation of structural stability. In 
particular they proved that linear flows with the same monodromy and winding can 
be deformed into each other. 



The tools forged above will now be used to compute the periodic Morse index in 
two different ways. The first is meant to illustrate the relation of the index with the 
stability properties of the flow, the second to recover the result of section 3.2.2. 



3.6. 1 The Conley and Zehnder index 

The most direct way to evaluate the periodic index is to consider intersections of the 
homotopy loop with the periodic pair {\2d , hd)- The restriction to regular crossings 
imposes some conditions on the linear flows. In particular flows with parabolic 
normal forms must be regularised by a positive definite perturbation in order to 
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render non singular the crossing form. 

With this proviso, the index is equal to the Maslov index of the third branch of the 
homotopy loop (230). This latter is most conveniently computed at s equal one for 
the family of symplectic curves (237): 

2d, hd), 

[0,T/2]} 

+ K{e(2*-^)J^fl,(l2,,l2d),[T/2,T]} (239) 

The first Maslov index is independent of the Lagrangian manifold associated to the 
periodic Bott pair (\2d , hd)- The statement is proven by the following reasoning. 
Since the matrix Pe(t) is periodic its index does not change if the interval [0, T/2] 
is replaced by [t', T/2 + 1'] for any t' such that no intersection occurs. Hormander's 
identity (223) and periodicity then show 



^Pe{2t), {U \2d): [0,T/2]} = ^Pe{2t),% [t\T/2 + t']} (240) 

for any 21 having no intersections with the graph of the flow at path end-points. In 
particular one can chose 21 ~ (hd^ — hd) and observe that any symplectic matrix 
admits a unique polar representation 



, Sy(2t) := (Pe^Pe)5(2t) 
Pe(2t) = (SyOr)(2t), ^ ^ / \ ^ (241) 

Or{2t) := [(Pe^Pe)-5Pe](2t) 

with both Or and Sy in Sp{2d). The polar representation can be embedded into a 
family of curves analogous to (237). Since Sy is positive definite, crossing occurs 
only when 



det[l2d+Or(2t)] = 



(242) 



as t ranges in [0 , 7/2]. The Maslov index can be thus computed from the intersec- 
tion of the orthogonal matrix Or(2 1) with the anti-periodic Bott pair. 
An orthogonal and symplectic matrix admits in general the block representation 



Or 



X Y 
-Y X 

and therefore it is similar to 



{X + i\)\X + i\) = \a 



(243) 



1 
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Id -i\d 




X Y 


1 


Id i\d 




Id i Id 




-Y X 




Id -i Id 





X + iY 
X-iY 



(244) 
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The analysis of the crossing forms yields modulo four the equality 



b^{Pe(2f),(Ul2d),[0,r/2]} = 

argdet[X(r) +iY(r)] - argdet[X(0) + ^ Y(0)] ^ ^^^4^245) 

The right hand side is an integer by periodicity. It coincides with twice the Gel 'f and 
and Lidskii [71] winding number of the linear periodic Hamiltonian flow Pe{2t). 

The only contribution to the second Maslov index in (239) comes when t is equal 
T/2. The crossing form can be straightforwardly evaluated from the normal forms 
of the stability blocks of the exponential flow. Following the same conventions of 
appendix 3.3 the Hamiltonian matrices associated to the blocks are 



HI ell. — 



U! 









H|hyp — 











HI 



lox. 



C^(l2 + J2) 

C^(l2 + J2)^ 



(246) 



A consequence of the analysis at the end of section E.2.2 is that inverse hyperbolic 
blocks behave at the origin as elliptic ones [59,100]. The simplest parabolic blocks 
have the normal form 



1 

In appendix (E.2.3) it is shown that (247) bifurcates under a generic positive defi- 
nite perturbation to an elliptic or to an hyperbolic block if, within the conventions 
adopted in this work, the sign of k is respectively negative or positive. In order to 
infer the sign of the infinitesimal characteristic frequencies of the elliptic blocks 
generated by the regularisation, one can reason as follows. A sufficiently small 
positive definite perturbation of the Lpcr. ([0, T]) does not affect the periodic Morse 
index because it removes the zero modes while preventing any eigenvalue from 
turning negative. Therefore the value of the Maslov index associated to parabolic 
blocks by such a perturbation of the Lper.([0, T]) must coincide with the minimum 
over the set of the Maslov indices generated by arbitrary perturbations of the same 
system [107,3,4,108]. This means that, once the Gel'fand-Lidskii winding number 
has been determined from the unperturbed system, the sign of the characteristic 
frequencies of elliptic blocks created by a positive definite perturbation must give 
rise to crossing forms with negative signature. A positive signature would produce 



F(^,T')L = 
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a larger value of the Maslov index. A rigorous proof of the above perturbative ar- 
gument can be found in [107]. 

Gathering the contributions of all blocks the resulting Maslov index is 

^^{e(^*-^)^^ ^ (b., [T/2, T]} = neu + n,^. - £ i^^^ (248) 

i=l 

where rieii, rii.h. and Up are respectively the number of elliptic and parabolic blocks 
in the monodromy matrix M = F(T, 0). 

The overall periodic Morse index is 

Up -. 

ind-Lpe,.([0, T]) = rieii + rii.h. - ^ + 2 tt) (249) 

i=i ^ 

with tn the Gel'fand-Lidskii winding number. Elliptic and inverse hyperbolic blocks 
can therefore be associated with the occurrence of an odd number of negative eigen- 
values in the spectrum of a periodic Sturm-Liouville operator Lper.. 

In the mathematical literature [136,137,142,1 10] the periodic Morse index expressed 
in the form (249) is often referred to as the Conley and Zehnder index. Conley and 
Zehnder [35] showed that the periodic Morse index (249) describes the foliation 
of the symplectic group Sp{2d) in leaves of different winding. Each leaf further 
decomposes in two connected components characterised by opposite sign of the 
determinant 



sign det[l2d - M] = (-ij'^-^eii-^ih. (250) 

In their paper [35], Conley and Zehnder analysed the stability properties of periodic 
orbits of non-autonomous Hamiltonian systems without parabolic blocks and in 
dimension larger than two. General proofs of the equivalence of the periodic Morse 
index with the Conley and Zehnder index were later given in [157,3,4] where the 
case of conservative systems with parabolic blocks is also encompassed. 

It is straightforward to infer the behaviour of the index under iteration of a prime 
periodic orbit. 

The winding number increases linearly with the number of repetitions r of the 
orbit. The sign (250) may instead change according to a nonlinear law. Namely 
under iteration elliptic blocks of frequency uj experience further intersections if 

r UJ becomes a multiple of 2 it. Inverse hyperbolic blocks generate a crossing form 
at each iteration since they are turned by even powers of the monodromy matrix 
into direct hyperbolic ones. Parabolic blocks do not change signature. Hence one 
concludes 
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ind- Lper.([0,rT]) 



l + 2int 




+ r fill,, - 



1 - sign Ki 



+ 2rro 



(251) 



i=l 



2 



The result holds for Hamiltonian systems without marginal degenerations. In [150], 
Sugita derived formulae (249), (251) by computing quadratic phase space trace path 
integrals on a time lattice. Sugita also pointed out that the negative sign in front of 
the parabolic block contribution renders in general the phase of the Gutzwiller trace 
formula independent from the signature of the energy block. A careful analysis of 
the behavior of the phase of the Gutzwiller trace formula predicted by (251) has 
been carried out in a recent paper by Pletyukhov and Brack [128] where several 
analytical and numerical examples are also investigated. 

More general explicit expressions of the periodic Morse index in the Conley and 
Zehnder form, encompassing marginal degenerations can be found in the paper 



The behaviour of Morse indices of Sturm-Liouville operators of variational origin 
was investigated in the topological setting in ref. [40]. The same results can also be 
obtained by means of pure algebraic methods [94,59]. 

Finally it is worth noting that the periodic Morse index can be identified with the 
winding number of the linear flow by applying the polar decomposition directly to 



3.6.2 Focal point description 

The periodic Morse index is needed in the Gutzwiller trace formula. Due to energy 
conservation at least one parabolic block is to be expected. Therefore in applica- 
tions [140] it is useful to exploit Hormander's identity (223) to derive the index 
without need of perturbative arguments. Inserting 21 ~ (l,^ © , © 1,^) in Duister- 
maat's formula (235) yields: 



The matrix D is specified by the block representation of the flow (166). Some te- 
dious algebra gives the expression of the concavity form: 



[108]. 



F(t,0) [105,140]. 



ind-Lper.([0,r])= Y nulD(i,0) 



0<t <T 



+ ind-Q(GrF(T, 0), I, © , © 1^; (bd, hd) ) - d (252) 



Q(GrF(r,0) ,(ld©0,0©l,);(l2d,l2d)) 

(CA-i)(r,o) At-i(r,o)-i, 
A-i(r,o)-id -(A-iB)(r,o) 



(253) 
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It is straightforward to see that in any vector basis where B(T, 0) and C(T, 0) vanish 
the index of (253) is equal to d. Thus the periodic index coincides with the zeroes 
encountered in the same vector basis by D(t, 0) through all the time interval [0, T]. 
More generally if B(T, 0) is nonsingular then one finds 

ind-Q(Gr F(T, 0), (I, © , © I,); {hd, hd) ) = 
ind-(D B"^ + B-^A - B"^ - B-^^){T, 0) + ind-[-(A-^B)(r, 0)] (254) 

Under the same hypothesis the use of Hormander's identity gives 
nulD(i, 0) = 

0<t<T 

^3{F, (0 © I, , © \d)} + d- ind-[-(A-iB)(T, 0)] (255) 

Since 

^3{F,(0©lrf,0©y} = ^nulB(t,0) (256) 

0<t <T 

the expression of the periodic index given in section 3.2.2 is finally recovered. 

The focal point representation can also be applied in a coordinate independent way. 
In particular it is possible to derive equation (249) (the Conley and Zehnder index), 
without need of the perturbative argument. 

Consider for example a periodic Sturm-Liouville operator Lper.([0, T]) such that 
the associated linear flow consists of a parabolic and an elliptic block. The Floquet 
representation is then 

F(i, 0) = Pe{t) e*-" ^ = Pepa..(i) e*-" © Peeii.(t) e*-" (257) 

As assumed above, phase cot of the exponential matrix of the elliptic block is sup- 
posed to range in [0, 2 7r[. 

The product and homotopy properties of Maslov indices illustrated in section 3.4 
permit to compute the periodic Morse index of Lper.([0,T]) as the sum of the 
Maslov indices of the stability blocks. Moreover, Hormander equality (223) can 
be separately applied to the periodic and exponential matrices entering the Floquet 
representation. The Maslov index of the periodic matrix is already known to coin- 
cide with twice the Gel'fand and Lidskii winding number vo. 
In order to compute the Maslov index of e*"'^ '^p''" , e*-*^ '^<=" in [0, T] one can apply 
Duistermaat formula (235) with 2t~(l©0,0®l). Consider first the elliptic block: 
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cos(a;t) 
- sin (cut) 



sin (cut) 
cos(a;t) 



(258) 



Crossings occur when 



det[(0 © 1) - 



;i ©0)] = -cos(tut) = 



(259) 



The condition is satisfied if cot is equal to 7r/2 or to 3 it/ 2. The corresponding 
crossing forms are equal to one. By equation (253) the order of concavity is 



ind Q(Gre 



Tjt H 



ind 



sin(a;T) 
cos(a;T) 

' -1 



cos(wT) 



0,0©l);(l2,l2) 
1 1 

cos(a;r) 

sin(a;r) 
cos{wT) 



1 = 



(260) 



It is equal to one for uT taking values in [0 ,it/2], is zero for uT in ]7r/2 , 37r/2] 
and finally minus one for uT in ]37r/2,2 7r]. Adding up the contributions of the 
crossing and of the concavity forms one has: 



^{gt Jt Heii.^ (1 © , e 1), [0, T]} = 1 



(261) 



independently of the value tuT in [0, 2 7r[. 

Turning to the parabolic block eq. (247) , it is seen that the only non zero contribu- 
tion comes from the order of concavity: 



ind-Q(Gre^-'^ ^^p-, (1 © , © 1); (b, I2) ) - 1 



ind 





— K 



1 — signre 



(262) 



Gathering all terms one gets into 



ind-Lpe,.([0, T[) = 1 - L_|g^ + 2ro (263) 

which is formula (249) applied to the example considered here. The block by block 
analysis of the Maslov index, has been used by An and Long [4] to prove the equiv- 
alence of the Duistermaat formula (235) with the Conley and Zehnder index. 
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3. 7 Phase space path integrals and infinite dimensional Morse theory 



Classical mechanics has its natural formulation in phase space. The lift of configu- 
ration space quantities to phase space was reiteratively used in the above presenta- 
tion of Morse index theory. It is therefore natural to wonder whether the treatment 
of the semiclassical approximation can be simplified by starting from phase space 
path integrals. 

Formally, phase space path integrals can be written directly in the continuum by 
replacing the quadratic kinetic energy term in the configuration space Lagrangian 
with its Fourier transform [51,101,135,165]. The propagator is then represented in 
the form 

K(Q, T\Q', T') = /'^^^"'^ V[q{t)pit)y^T' dtWr-H{p,q)] (264) 

Jq{T>)=Q' 

The domain of integration is asymmetric between position and momentum vari- 
ables as the boundary conditions impose a constraint only on position variables. 
More general phase space path integral expression are obtained if (264) is used to 
compute the time evolution of quantum observables other than the propagator. 
Formally the derivation of the semiclassical approximation of section 2.2 goes 
through also for phase space path integrals. In the presence of boundary conditions 
03 of the form (58) it leads to the quadratic action 

S^S = J^dtSx'^ZafsSx'^ , a,P^l,...,2d (265) 

where 6x = {5q, 5p) denotes the fluctuations around the phase space classical tra- 
jectory Xci while 



Zap '■— Ja/3 ^ — (266) 

with Hajs the Hamiltonian matrix (147) at r equal zero. The operator 3<b is then 
self-adjoint with respect to the phase space scalar product 

{x, y) = Ij dtx-5^,py^ = £ dtx^y (267) 

The scalar product presupposes a suitable rescaling of positions and momenta. The 
fact is not at variance with classical mechanics where position and momenta lose 
their physical interpretation under canonical transformations [7,98]. 
Forman theorem can be applied to recover the configuration space result 
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I Deta<8( [T', T] ) I = I x<8 det [Yi + Y2 F(r, T')] \ (268) 

The operator is a first order Dirac operator [101,143,137]. At variance with the 
Sturm-Liouville operators of the Lagrangian formulation the eigenvalue spectrum 
is unbounded from below and from above. Hence the Morse index cannot be de- 
fined as in configuration space. However, the lattice approximation identifies the 
index of any non singular discrete version of 3<b as the difference between the pos- 
itive and negative eigenvalues. It is possible to prove [138,139] that in the infinite 
time lattice limit such quantity converges to the Morse index found above work- 
ing in configuration space (see also [110]). Difficulties may arise if one works in 
directly in the continuum. The discussion can be made more concrete for periodic 
boundary conditions. The eigenvalue spectrum of limit fluctuation operator 3per. is 
invariant under periodic transformations in Sp{2d) continuously connected to the 
identity: 



S^aS - J ^ - (S^ HS - JS) (269) 
at 



The choice 



S(t,r') = F(t,r')e-"^(*-^') (270) 

with H defined as in (236) replaces the Hamiltonian matrix in Zpct. with its aver- 
age over one period. The eigenvalue problem is then explicitly solvable in general 
[44,45,54,150]. The result is evidently independent of the winding number of the 
monodromy of the linearised flow. This means that the explicit diagonalisation of 
Zpci. will give only the correct sign of the functional determinant but not the Morse 
index. The semiclassical approximation needs the square root of the determinant. 
Thus, in the phase space continuum limit the phase factor of the path integral can- 
not he defined by simply diagonalising the fluctuation operator. Instead, it is defined 
by the the index of the eigenvalue (Fredholm) flow of the homotopy transformation 
between and another Dirac operator whereof the index is a priori known. This 
is the content of the infinite dimensional Morse theory developed in ref.'s [143,137] 
where the topological methods of sections 3.4, 3.5 are systematically applied. Thus, 
the evaluation of the periodic Morse index turns out to be an example in the sim- 
plest physical context of a general Fredholm flow theory which has found in recent 
years wide applications in field theory starting with the treatment of the SU{2) 
anomaly given by Witten in [160]. Recently Witten's approach has been applied by 
Sugita in [150] to the derivation of the phase of the Gutzwiller trace formula. 

Once the more subtle continuum limit rules of calculus are understood, phase space 
path integrals offer the advantage of an elegant formalism able to deal with a wider 
class of problems than their configuration space homologues. A typical case is the 
quantisation of systems classically governed by singular Lagrangians, see for ex- 
ample [60]. 
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4 Stationary phase approximation of trace path integrals 



If the Lagrangian admits periodic extremals and is invariant under the action of a 
Lie group, the periodic Sturm-Liouville operator associated to the second variation 
is singular. The quadratic approximation of chapter 2 diverges. As it is well known 
in soliton and instanton calculus [17,130,131,132,152,153,34,154,55], the diver- 
gence can be cured using the Faddeev and Popov method [61,60]. Gutzwiller's 
trace formula follows from a direct application of the method both for Abelian 
or non Abelian Lie symmetries. The expressions of the degenerate trace formu- 
lae derived by Creagh and Littlejohn [36,38] with traditional WKB techniques are 
recovered in the more general path integral formalism. 



4. 1 Trace path integrals and zero modes 



In chapter 2 the semiclassical approximation of path integrals was derived under 
the hypothesis of a non singular fluctuation operator. If the action is invariant under 
time translations, the trace path integral 



TrK= I V[^q{t)]e^^ (271) 

violates the assumption. A time translation maps a solution of the classical equa- 
tions of the motion into another solution. Thus, the set of trajectories on a periodic 
orbit defines a degenerate extremal of the action. The infinitesimal generator of time 
translations is the time derivative d/dt. The induced vector field along a T-periodic 
trajectory qci{t) is the velocity qci{t) which is therefore also a periodic Jacobi field 
for the second variation. To see this, it is enough to differentiate the Euler-Lagrange 
equations specifying the periodic trajectory qce{t): 



d 
di 



dC 



d dC 



dq'^{t) dtdq'^{t) 







(272) 



and observe that the result can be rewritten as 



(Lper.gc£)W = (273) 

Hence the velocity field is a zero mode of the second variation operator Lp^r. around 
the periodic extremal qci{t). 

Nother theorem provides a general mechanism to associate zero modes to con- 
tinuous symmetries of the action. Let for example G be an unbroken non-trivial. 
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compact, connected finite-dimensional group G of symmetries of the action con- 
stituted by isometrics which also leave invariant vector and scalar potentials. Fur- 
thermore, assume G to be globally parametrised by coordinates (t^, t^). In such 
coordinates the origin is supposed to coincide with the identity transformation. Any 
dynamic extremal qci{t) of the path integral action which breaks completely G is 
embedded in an iV-parameter family of solutions of the Euler-Lagrange equations. 
The generic element of the family is obtained from qct{t) by the symmetry trans- 
formation 



?S[t^...,t-](^):=r(?d(^),t^...,0 (274) 

The Jacobi fields are spawned by the vector fields induced along the trajectory by 
the generators of infinitesimal transformations around the identity: 



m = lim V°feW.O.....Q.^^.O....,Q)-tf<W (275) 

The Jacobi fields are T-periodic if qci{t) is T-periodic. Thus they are zero modes 
of Lper., the fluctuation operator of qci{t)- The velocity field can be as well encom- 
passed by the definition (275). This is done by shifting t to t + 1° and by re-defining 
the transformation law </? as the composition of the group action on the trajectory 
coordinates with the flow solving the Euler-Lagrange equations. The set of param- 
eters (t*^, t^), collectively denoted by t, are the moduli or collective coordinates 
of the degenerate stationary point. 

The occurrence of zero modes among the quadratic fluctuations is not a path in- 
tegral peculiarity. The same problem is present already for ordinary integrals. A 
classical example [165] is the "zero-dimensional" field theory 



i{h) = /' df^Qe-'-^ (276) 



with action 



U{\Q\)^U{\RQ\) {111) 

invariant under the group of (i-dimensional rotations SO{d) and growing to infinity 
for large value of the radial coordinate 

In the small h limit the integral can be performed by means of the steepest descent 
method. The first two derivatives of the exponent in Cartesian coordinates 
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aQ"^'^'^ \Q\d\Q\ 

ag^ag/^^'"^'^" igi^ d^igi^''^'^^ |g| [''^^ |gp 



(|g|) (278) 



d\Q\ 



show that quadratic fluctuations around the any extremum Q* different from zero 
are governed by the projector 



dQ^dQf^'^^ " |g*|2 d?\Q\ 

The projector has d — 1 zero eigenvalues. The number of zero modes is equal to 
the dimension of the quotient SO{d)/SO{d — 1) between the original group of 
symmetry and the subgroup 5*0(0? — 1) of transformations leaving an extremum 
invariant. The phenomenon is general and well known in field theory as it provides 
the mechanism for spontaneous symmetry breaking. The subgroup Gs of a symme- 
try group G leaving invariant a "vacuum state" is then called the stationary group 
or stabilizer of the vacuum. The dimension of the quotient V = G/G^, yields the 
number of Goldstone modes, emerging from the symmetry breaking [77,78,165]. 
The breakdown of the quadratic approximation is due to the violation of the orig- 
inal symmetry of the action by its extremal point. The problem is obviated by the 
use of the symmetry to eliminate the ignorable degrees of freedom from the ap- 
proximation through a suitable change of variables. In the example (276) this is 
simply done by turning to spherical coordinates. The steepest descent - stationary 
phase approximation is then performed on the reduced action while the ignorable 
degrees of freedom are integrated non-perturbatively. 

The same method can be applied to path integrals. For example in ref. [133] the 
Edwards and Gulyaev construction (appendix D) was used to write quantum me- 
chanical and field theoretic trace path integrals in spherical coordinates. The en- 
ergy density of systems governed by Lagrangians invariant under rotations was 
then evaluated by applying the semiclassical approximation to the radial coordi- 
nate. The original symmetry is thus always restored in the final expression of the 
steepest descent approximation for all systems with a finite number of degrees of 
freedom. In other words, no spontaneous symmetry breaking can occur in quantum 
mechanics [126]. 



4.2 Faddeev and Popov method and moduli space 



Performing as in ref. [133] the removal of ignorable coordinates by change of vari- 
ables in the path integral may turn to be unwieldy. A viable alternative is provided 
by the Faddeev and Popov method [61] originally devised for the quantisation of 
gauge field theories. In order to streamline the discussion, the hypothesis is made 
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that the dynamic extremals of the path-integral break completely the G x [0, T] 
symmetry of the trace path integral. 

Symmetry transformations acting on a path in the loop space define an equiv- 
alence class of configurations or group orbit which leaves the action functional 
invariant. A periodic degenerate stationary point of the action represents an equiv- 
alence class of periodic classical trajectories. Hence the path integral is expected to 
be proportional to the volume of the group orbit. For a compact group G x [0, T] 
such volume is finite and it is natural to extract its contribution before proceeding 
to the semiclassical approximation. In order to do so, one recipe is to restrict the 
path integral to an "hyper-surface" in which intersects any group orbit only 
once. This means that if 



Ta[q{t)\^Q, Va = 0,...,7V (280) 

is the equation specifying the hyper-surface, then the equation 

^a[V'(gW,t)]=0, Va = 0,...,7V (281) 

has only one solution in G x [0,T]. The path integral constrained to the hyper- 
surface ranges then over LSDt/ (G x [0, T]) and does not contain any longer the zero 

modes produced by the symmetry. 

The decomposition of any closed path in [0, T] into a classical T-periodic trajectory 
and quantum fluctuation is accomplished by setting 

?"W=?c^[t]W + ^<^?wW (282) 

The quantum path on the left hand side of (282) is by definition independent of the 
moduli t. The simplest choice of the constraint T is then 



:FaW)] [ dt[q-{t) - q'^^,^,^{t)f-^{t) , Va = 0, 7V(283) 

Applying Lie's first fundamental theorem (appendix F) the derivatives of the classi- 
cal trajectory with respect to the moduli are expressed in terms of the Jacobi fields 
(275) induced by the infinitesimal generators of the group transformations: 

'%^{t)-<{i)'§f{t)i[t) (284) 

The matrix 9^5 (1) depends only on the moduli and coincides with the identity when 
t is equal to zero. The constraint (283) is inserted in the path integral in the guise 
of the Faddeev-Popov expression of the unity 
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/ TT de det I 

iGx[0,T]f4 t 



di'> 



N 



j=0 



(285) 



The matrix elements of the Jacobian, 



^ = /;* {rw - fi„(t)i%^w - ^w%^w) (286) 

are computed on the intersection of the group orbit with the manifold (280). Once 

(285) has been inserted in the path integral the order of integration between the 
moduli and the path integral measure is exchanged. Quantum fluctuations are then 
written as the -series 



<^?MW = EcnAifl"(t) (287) 

n 

extended over the eigenstates of the fluctuation operator Lp}j_.^ around the classi- 
cal trajectory q"i[q{t). The geometrical meaning of the constrain (283) becomes 
manifest: it restricts the span of quantum fluctuations to "massive" eigenstates, or- 
thogonal with respect to the L^-scalar product to the nuUspace of Lp}^. . Actually, 
an orthonormal basis for this latter is composed by vector fields of the form 



Ail%W = {Z-'t -j^it) £{t) , a = 0, TV (288) 

dq'ce 

as it follows from the use of the Gram-Schimdt recursion method. Since the sym- 
metry is an isometry for the metric. 



9a^i<l)-g^Q^9.sm,i)) (289) 

it is enough to compute the triangular matrix around the identity transformation 
using the Gram-Schmidt recursion relations: 



K.m.,l{t) 

etc. 



m 

Oo , Jo) 



Oo Jl)^ 
J0> 



(290) 
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The path integral measure was argued in chapter 2 to be proportional to the infinite 
product of integrals over the amplitudes c„ of the L^-expansion of the quantum 
paths. Therefore, carried under path integral sign (285) becomes 



N 



Gx[0,T]f4 



det 



N 



3=0 

N 



. N N 

JGx[0,T]^^Q • Q 



(291) 



with a,b — 0, N. The 5-functions remove the divergences from the semiclassical 
approximation: 



TrK = 

y /-^ f di'^dG [ VCUdcJ det Z(°) TT S(ca) + o(Vh) (292) 

The sum over o ranges over all the distinct degenerate stationary points of the ac- 
tion. The measure dG 



N 



dG= lldtadet9i{t\...,i^) 



a=l 



is invariant over the compact group G ([19] and appendix F). 
The path integral (292) yields 



(293) 



TrK 



5(0) 

oSp.o. 



7 dt^dG 

JGx[Q,T] 



md-4°' ([o,T])+^ 



Deti4°>J[0,T]) 



det Z(°) 2 



(294) 



Det_L is the functional determinant produced by the restriction to massive modes. 
Symmetry transformations are by (289) isometrics also for the L^-scalar product. 
The integrand in (294) is therefore function only of the degenerate stationary point 
of the action. One can rewrite 



TrK 



E 

oep.o. 



T\G\ 



' ft 



ind-4°l([0,T])+^ 



Det L' 



.([0,T]) 



det Z(°) 2 



(295) 



with 
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dG 



(296) 



The factor [2 71 he'^)^~ in the denominator of (295) comes from the definition 
of the path integral measure (compare with appendix H). Namely the integrals 
over the amplitudes of the expansion of the quantum paths are normalised to 
(27r/ie^i)i 

In the following section it will be shown how the use of Forman's theorem recasts 
the last formula in terms of invariant quantities of the classical periodic orbits. 

An historical remark before closing the section. The Faddeev-Popov method was 
applied to the treatment of constrained quantum mechanical system by Faddeev 
himself in ref. [60] using phase-space path integrals. However, in ref. [131], Poly akov 
credits Bogoliubov and Tyablikov [27] to have been the first ones to treat non Gaus- 
sian fluctuations emerging in the steepest descent - stationary phase approximation 
with the method expounded above. 

The Faddeev-Popov method has found wide application for the quantisation of 
quantum, statistical mechanical and field theories in a soliton or instanton back- 
ground. Classical references are [73,75,33,17,130,131,152,153,34,154]. A recent 
review is [55]. The treatment of collective coordinates in problems with dynamic 
extremals only partially breaking the symmetry group has been analysed in [74]. 
The application to scattering problems is discussed in [156]. 
The consistency of the Faddeev-Popov method with a systematic asymptotic per- 
turbation theory beyond leading order has been checked for example in references 
[2,161,6]. The investigation of higher order corrections is sometimes simplified if 
the Faddeev-Popov method is implemented by means of ghost fields and the impo- 
sition of BRST conditions [16]. The procedure is reviewed in [25,70,165]. 

4.3 Zero mode subtraction 

The explicit expression of the functional determinant of the massive modes can be 
computed adopting the same strategy pursued in the evaluation of the Morse in- 
dex of Sturm-Liouville operators with zero modes. A positive definite T-periodic 
generic infinitesimal perturbation rU is introduced to subtract the zero modes. For- 
man's theorem can be applied directly to the non-singular operator Lper. ,t([0, T]). 
Finally, the functional determinant of massive modes is recovered by taking the 
limit 



The product in the denominator runs over the eigenvalues of Lpcr. ,r ( [0 , T] ) flowing 
into zero modes at r equal zero. In order to evaluate (297) it is enough to compute 
the eigenvalue product up to first order in r: 



Det^Lpe,.([0,r])| =lim 



det[l2d-F,(T,0)] 



(297) 
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fi 4,. = r^-^^ ^fls^r^^lx + o{T^^') (298) 

The denominator cancels out with the Gram-Schmidt Jacobian det Z. The cancella- 
tion is not accidental since the functional determinant restricted to massive modes 
by definition cannot depend on zero modes. For the same reason the averages of 
the perturbation over the periodic Jacobi fields should cancel out. The mechanism 
of the cancellation is governed by the parabolic blocks of the monodromy matrix 
of the periodic orbit. The monodromy matrix is obtained from the linearisation of 
the classical flow around any trajectory q'^^{t) on the periodic orbit by setting 

M = F(r, 0) (299) 
The monodromy matrix admits an orthogonal decomposition 



M = M" © (300) 

with M", M-*- governing respectively the stability of degrees of freedom longitudi- 
nal or transversal to the orbit in phase space. The decomposition corresponds to a 
partial reduction to normal form of the monodromy. The continuity of the eigen- 
value flow of a symplectic matrix under continuous parametric perturbation permits 
to write at r different from zero 



= M| ® (301) 

The perturbation shifts the stability of the blocks away from the marginal case. 
Hence, the block M l[ pairs up to degrees of freedom longitudinal to the orbit only 
at T equal zero. The functional determinant factorises as 

detihd - M,] = det[l, - Mil] detihd-i - M^] (302) 
with / < 2 (i the dimension of the parabolic subspace of phase space. The limit 



limdet[l2d_^ - M^] = det[\2d-i - M^] := det^i^d - M] (303) 

T— »0 



is non zero and independent from the moduli. Linearised flows along different tra- 
jectories belonging to the same degenerate extremum are related by similarity trans- 
formations which leave detj_[l2d — M] invariant. 



4.3.1 Abelian symmetries 

The symmetry group is Abelian if the generating functions of the symmetry trans- 
formations are in involution (appendix E.1.5): 
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{Ha, "Wblp.b. = , 



< a,b < N < d 



(304) 



The Poisson brackets (304) are equivalent to skew-orthogonality relations between 
the phase space lifts 



Ja{t) = 



(Vj„)(t) 



(305) 



of the N + 1 periodic Jacobi fields. By assumption the Jacobi fields are linearly 
independent. Hence + 1 is at most equal to d. The construction of a symplec- 
tic basis requires the introduction of other + 1 linear independent generalised 
eigenvectors of the monodromy matrix 



M 



fulfilling the normalisation conditions 



Q,..,N 



(306) 



a,b = Q,..,N 



(307) 



The generalised eigenvectors are mapped by the symplectic matrix J into a basis 
dual to the Ja's. By means of a linear symplectic transformation it is also possible 
to reduce M" to the orthogonal product of + 1 two-dimensional parabolic blocks 
of the form given in appendix E.2.2. Otherwise phrased, involution renders the 
normalisation conditions (307) compatible with the existence of a reference frame 
where each of the configuration space projections of the periodic eigenvectors of 
M coincides with one of the iJ'-orthonormal zero modes of Lpcr.. Introducing an 
infinitesimal periodic perturbation, the same calculation of appendix E.2.3 yields 



M 



\n+i V 



(308) 



The terms of order 0{t^^^) or higher neglected in (308) do not contribute to the 
limit 



lim 



det[l2(Ar+i) - Ml 



r^+i det{(j,,Uj,)} 



detVl 



(309) 



The determinant on the right hand side is just the product of the non-diagonal ele- 
ments of the normal forms of the parabolic blocks of the monodromy matrix repre- 
sented in the canonical symplectic basis 
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TV 

|detV| = I n^^il (310) 

i=0 

Physically 's are the characteristic growth rates of the linear independent marginal 
instabilities of the degenerate periodic orbit. The functional determinant is therefore 
independent on the moduli. Provided all the orbits contributing to the stationary 
phase are prime, the semiclassical trace path integral with Abelian symmetries is 
within leading order 



^1^1 ^Y-¥(-d-41([0,T])+^) 

TrK^ J2 M+r , = (311) 

oep.p.o. {2 7rh)— J\ det Vo det±[l2d - MJ | 



Were the stationary phase on a time interval r T with r integer dominated only by 
the iteration of prime orbits, the above formula would admit the extension: 



^ (ind-41([0,rT])+^) 

TrX^ V V , = (312) 

oep.p.o. r-=i {27rhr) — J\ det Vo det±[\2d - K]\ 



smce 



M(rT) = M'^(T) (313) 
4.3.2 Non Abelian symmetries 

A symmetry group is non- Abelian if the generating functions {Tiaj^i of group 
transformations form a non trivial Poisson brackets algebra (appendix E) 

{^a,^6}p.b. = -C^.Hc + D,,, a, 6=1, ...,7V < 2d (314) 

In consequence, the skew products among the right periodic eigenvectors of the 
monodromy matrix 



Jr = J^"''^, a = 0,...,N (315) 

with Ho coinciding with the Hamiltonian H, are now fixed by the Lie algebra of 
the group: 



Ja -^aP Jb = -V^a , '^bjp.b. (316) 
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For any trajectory on a generic orbit of period T, the zero eigenvalues of the result- 
ing + 1 X + 1 antisymmetric matrix correspond to the existence of Casimir 
operators in G x [0, T]. The Hamiltonian H furnishes an obvious example of a gen- 
erating function commuting with all the elements of the Poisson algebra. 
Let A; -I- 1 be the number of eigenvectors skew orthogonal mutually and to the re- 
maining N — k. Asin the Abelian case the diagonalisation of the block monodromy 
matrix spanned by such eigendirections requires the introduction of k + 1 skew- 
orthonormal generalised eigenvectors 



^ " " W a ^ a = 0,...,k , b = 0,...,N (317) 

The remaining N — k eigenvectors have non-degenerate skew products. From their 
linear combinations it is possible to form a symplectic basis for the linear subspace 
they span. The limit 



lim 

TlO 



det[l^+fc+2 - 



T^+i det{0,,Uj,)} 



detWdetVl (318) 



yields two contributions. The determinant det V retains the same definition as in the 
Abelian case. The new term det W accounts for the different normalisation, dictated 
by (316), of the configuration space projections of the periodic eigenvectors of the 
monodromy matrix from the corresponding Jacobi fields describing the zero modes 
of Lper. . In consequence, det W coincides with absolute value of the inverse product 
of the non-zero eigenvalues of the antisymmetric matrix specified by (316). If only 
prime periodic orbits contribute, the leading order result is 



r \r\ /-^-"^ (ind-41([o,T])+i^) 
oep.p.o. (27rn) — J|detWodetVodetx[l2d- MJ| 



The result is illustrated by an example common in applications. 

Consider an action functional with time autonomous Lagrangian invariant under 
SO{?>). For the sake of simplicity in what follows the Darboux variable x — {q,p) 
describes a six dimensional phase space. The hypothesis is non restrictive as it cor- 
responds to the choice of a frame of coordinates such that x describes the parabolic 
degrees of freedom in the monodromy matrix. The first integrals associated to the 
symmetry are the Hamiltonian and the three components Ala> a = 1, 3 of the 
angular momentum with Poisson algebra 



{Ma,Mb]^.^.=el^,Mc, a,6,c= 1,...,3 (320) 
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with Eabc the completely antisymmetric tensor. Let Ra, (a = 1, .... 3) denote a rota- 
tion around the a-th Cartesian axis. A canonical parametrisation [144] of S0{3) is 
provided by the angular momentum versor 



'^■■;^-|l^-'' = R.(t')R3(a (0,0.1)' 

:= (cos sin t^, sin sin cos t^)^ (321) 

and by the angle of rotation around the direction of the angular momentum. In 
terms of the moduli (t^, t^, t^), the generic element of SO (3) admits the represen- 
tation 

R{t\e,e) = R3(t^)R2(t')R3(t')[R3(t')R2(t')]^ (322) 

with 

0<tSt^<27r, < < n (323) 

By (284) the Jacobian entering the definition of the invariant measure can be ex- 
tracted from the left invariant differential over SO (3) 



[R-'dR){t\ e, t^) = (it"<(t\ t^ 



(324) 



The differential is written in matrix notation in the basis of the infinitesimal gener- 
ators of rotations around the three coordinate axes 



dRa 

ds 



(325) 



s=0 



The invariant measure reads 



dS0{3) = n di^"^ det 9l(t\ t^ t=^) = [] di^"^ 4 sin t^^) sin^ — (326) 

a=l a=l ^ 

The result does not depend on the use of the left-invariant differential (324) to con- 
struct the measure since SO{3) is compact and connected: the invariant measure is 
therefore unique [76,144]. 

The squared modulo of the angular momentum defines the generating function as- 
sociated to the Casimir operator of 5*0(3). It is therefore convenient to choose the 
periodic eigenvectors of the monodromy matrix in the guise 
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J^-^'-'n^^-P^ (327) 



with a, 6 = 1, 3. The {na}^^^ are time independent orthonormal versors in 
such that 



s: - 6' r 



(328) 



The explicit form of the versors, modulo a normalisation factor, can be derived 
directly from the expression (284) of Lie's first fundamental theorem applied to 
SO{?>). With the choice (327) the last three periodic eigenvectors have projections 
in configuration space mutually orthogonal with 77,3 always pointing in the direction 
of the angular momentum. 

The Poisson algebra of the generating functions yields the skew orthogonality rela- 
tions 



J:UJS = --Mc^L^a^ft = -WMWhle^aeuinl (329) 

The diagonalisation of the monodromy matrix requires the introduction of two gen- 
eralised eigenvectors 



dT 

M-^JSM\\-J\\M\\ + ^\\M\\Jr (330) 



where h\\m\\ physically gives the variation of the total angle of rotation around the 
angular momentum versus a change of the absolute value of the angular momentum 
itself. Finally for any generic periodic orbit described by the system one can apply 
(319) with 



detW| 
detVI 



\\M\\' 
dT 



(331) 
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4.4 Gutzwiller trace formula 



Gutzwiller's trace formula is a direct application of the general theory of path inte- 
gration over loop spaces. 

The relation (16) between the energy spectrum of a time autonomous quantum 
mechanical system and the propagator is amenable to the path integral expression 



roo ZT r 

p(E) = - lim Im / -e'— / V[^q(t)]e^ 

The semiclassical approximation must be performed on the Fourier transform of 
the propagator trace. The integration over the time variable entails the evaluation 
of the propagator for infinitesimally small times T. At variance with the finite T 
case, in such a limit the semiclassical approximation cannot be identified with the 
stationary phase approximation. The propagator becomes proportional to a 5-Dirac 
distribution for zero time increments: 



(332) 

E=R.e Z 



K,{Q,T\Q',0) ~ (-^!^)%-A/o"'^*^^fe'^*) (333) 
\2t:z nT J 

The divergence of the prefactor does not allow to neglect it in comparison to the 

fast variation of the phase for h going to zero. Thus, the stationary phase cannot 
be directly applied in this limit. In consequence, the semiclassical energy spectrum 
consists of the separate contribution coming from short orbits with period tending 
to zero and from long orbits of finite period T. 



4.4.1 Short orbits 

The short orbit contribution to the trace formula was investigated in detail long ago 
in the review by Berry and Mount [23]. 

The short orbit contribution is obtained starting from the Fourier transform of the 
short time representation of the propagator on closed paths 



The insertion of the short time approximation into the energy density (332) provides 
in an Euclidean space the asymptotic expression 
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r ^ f d'^Pd'^Q f°° Z-H(.P,Q) J, vT 

P3,oX^)~ lim Im— / — — — / ^ ^ 



imzio TrhJR^'i {27rhy Jo 

The short orbit contribution brings about a microcanonic average over the energy 
surface which dominates the asymptotic of the energy density at large values of E 
[21,126]. 

In the typical case of a kinetic plus potential energy Hamiltonian 



^ = E^ + ^(Q) (336) 
^ 2m 

the short orbit contribution reduces to 



the prefactor being the measure of the surface of a unit sphere in d-dimensions. 



4.4.2 Long orbits 

The long orbit contribution is just the Fourier transform of the semiclassical asymp- 
totics of trace path integrals. The time-energy Fourier transform makes sense be- 
cause generic conservative Hamiltonian systems are expected to exhibit periodic 
orbit in smooth families versus the energy. The time-energy Fourier transform in 
the semiclassical limit reduces then to a further stationary phase approximation. Al- 
ternatively, it is also possible to apply the stationary phase approximation directly 
to the energy density formula (332). In such a case the functional to extremise is 



W{q,q)^ET+ [ dtC{q,q) (338) 
Jo 

with the loop space condition 

q'^{t + T)^q'^{t) yt (339) 

For differentiable paths this latter entails 



^{t + T; T, ...) = -q^{t- T, ...) + ^{t; T, ...) (340) 
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the partial derivative with respect to T affecting only the parametric dependence of 
the path on its period. 

On a T periodic trajectory the functional (338) is equal to the classical reduced 
action [7,98]: 



>V(?c€,gc^) = j^ {dq^^m (341) 

The reduced action is function only of the orbit and not of individual trajectories 
wending their path on the orbit. By (339) and (340) the functional is stationary 
when 



Jo 

dC d dC 

DC— 

" ■ dq'^ dtdq'^ 



C-cf 



dC 
dq'^ 



t=T 



ol 



(342) 



vanish. Extrema are periodic orbit of now fixed energy E. The first of the equations 
(342) establishes therefore an energy-period relation E = E{T). 
Evaluated over any classical periodic trajectory q^, the second variation does not 
contain mixed terms since drqa is a (non-periodic) Jacobi field: 



dE 

S'yV{qAt),qceit)) = rdt5q^{Lp,,.U^5q^ 



(343) 



The integral over the fluctuations finally yields: 



Pl.o.(-S) ~ 

Im V V — — (344) 

h ^ ^ (2 1 1^ det W„ det V^^^ det^hd - I ^ 

where the sum ranges on all orbits "o" at energy E and their r-th iterates. The phase 
factor associated to extrema of the reduced action is 



K,r - ind-L£. [0, r TJ + 1(^1 + sign^ j + ^ (345) 

The use of the Conley and Zehnder representation of the periodic Morse index (25 1) 
shows that the signature of the energy block cancels out from the phase, [150]. If 
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the energy is the only conserved quantity, the formula further simplifies (see also 
appendix E.2.3) 



det W„ = 1 



IdE, 



r dT 



^detVW 



(346) 



so to retrieve Gutzwiller's original result given in chapter 1. The opposite limit of a 
separable system is treated by subtracting the d zero modes associated to the prime 
integrals in involution. In this way the result of Berry and Tabor [24] is recovered. 
The reader is referred to the existing literature [133,54,36,37] for a detailed discus- 
sion of the semiclassical quantisation of the extra integrals of the motion. 



4.4.3 Example 

The simplest application of the trace formula is the derivation of the semiclassical 
energy spectrum of a one dimensional particle with kinetic plus potential classical 
Hamiltonian. Moreover the potential is supposed to be a single well. 

The short orbit contribution is proportional to the period of the accessible orbits at 
energy E 



Qmax dQ 



m 



ToiE) 



Turning to long orbits one observes that (345) is equal to 



(347) 



K., = 2r 



(348) 



as it follows from the analysis carried out in sections 3.3 and 3.6. Thus the oscillat- 
ing term in the trace formula is 



oep.p.o. ^ r=l 



'rW{E) 

h 



7rr 



(349) 



Gathering the short and long orbit contributions gives 



P{E)= E 



oep.p.o. 



- E 



To{E) 

TT h 

UE) 



1 + 2E 



rWiE) 

cos I : TT r 

r=l \ 



E 



oep.p.o. r=—oo 



/VV%B)_1 
2wh 2 



(350) 
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The series over r yields a discrete representation of a train delta functions peaked 
at 



for n an arbitrary integer. Thus the Bohr-Sommerfeld [28,148] quantisation rule is 
recovered. 



4.4.4 Non generic degenerations 

Zero modes occur generically in correspondence to continuous symmetries of the 
action. However, it is not infrequent to encounter in applications marginal cases 
where unit eigenvalues of the monodromy matrix do not stem from the invariance of 
the action under a Lie group. The zero mode subtraction method expounded above 
proves nevertheless useful in order to compute the functional determinant of the 
massive modes. The contribution of the degrees of freedom associated to marginal 
zero modes can be then computed using approximations higher than quadratic. This 
is done, for example, by finding the normal form around the periodic orbit of the La- 
grangian and then retaining only the leading order projection along the zero mode 
eigendirections. The procedure is discussed in details by Schulman in his classical 
monograph on path integration [145]. 
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5 Conclusions 



The derivation of the Gutzwiller trace formula is notoriously difficult. 

In the author's opinion the difficulty is considerably mitigated by the use of path 

integral methods. The opinion is based on two strictly intertwined reasons. 

The first reason is that trace path integrals offer a global, canonically invariant ap- 
proach to the semiclassical approximation of the energy density of a quantum sys- 
tem. This is in contrast to the WKB methods usually applied in the literature of 
the Gutzwiller trace formula. Broadly speaking, WKB methods proceed by repeat- 
ing explicitly in the case of the trace of the propagator the construction made in 
section 2.1. The phase of single orbit contributions is determined by keeping track 
of the caustics encountered along a reference trajectory and finally computing the 
order of concavity when the trajectory has covered the entire orbit. 
Path integrals permit to concentrate directly on the properties of the propagator 
trace. The change of point of view is of particular advantage in the interpretation 
and then practical computation of the phase factors of the orbit contributions. Here 
the difference between path integral and WKB methods can be summarised as a 
shift of the focus from the local geometrical properties of classical phase space to 
the global topological properties of second variation self-adjoint operators. Morse's 
theory of variational calculus in the large provides all the information to compute 
the index associated to periodic orbits. It is worth stressing that the topologically 
invariant Morse index for closed extremals was already computed by elementary 
methods in Morse classical monograph [118] (see also [119]). 
The Lagrangian manifolds techniques later developed by Bott, Arnol'd and Duister- 
maat establish an extremely powerful connection between Morse index theory and 
Fredholm flow theory. Local geometry of phase space then re-emerges but only as 
a result of topological invariance. 

The index of the second variation defines as well the index of the functional deter- 
minant of the self-adjoint operator associated to the second variation. Whenever the 
classical kinetic energy is strictly positive definite, configuration and phase space 
path integrals provide equivalent representations of the propagator trace. In conse- 
quence the Morse index for closed extremals defined in configuration space must 
coincide with the Conley and Zehnder index [35,143,137] associated to the Dirac 
operators governing the second variation in phase space. This latter is widely used 
in recent mathematical investigation (see overview in [142] and the monograph 
[110]) aimed at establishing the general conditions which guarantee the existence 
of periodic orbits on a manifold equipped with a symplectic structure. 

The second reason dwells in the general use of the path integral methods wielded in 

the derivation. They provide a common language for the treatment of quantum and 
statistical finite dimensional and field theoretic models. In this framework, beside 
its intrinsic importance, the Gutzwiller trace formula acquires also a valuable ped- 
agogical significance. It provides a paradigm for the application of general math- 
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ematical ideas which have proven of large use and relevance in different parts of 
physics. The author's hope is that the present may also serve as an illustration of 
such ideas accessible to a broad physical and mathematical audience. 
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A Elementary geometric concepts 



A smooth, torsion-free manifold with strictly positive symmetric metric tensor gap 
is said to be a Riemann manifold. The metric with its inverse are used to lower and 
raise indices 



v^^g-^vp (A.1) 
The covariant derivative of a vector field 



- ^ (A.2) 

is compatible with the metric if for any pair of vector fields v"', x" evaluated along 
an arbitrary curve G Tl one has 



The identity is always satisfied if the connection Tg satisfies 



r^. = (A.4) 

On a Riemann manifold the above compatibility condition is uniquely solved by 
the Christoffel symbols : 



^lu = + du9p^. - dpg,,) (A.5) 

A non trivial metric tensor may arise from the parametrisation of an Euclidean 
space in non Cartesian variables. The curvature tensor, the commutator of two co- 
variant derivatives, discriminates between this latter case and that of genuinely non- 
Euclidean space 



^^^.^'':=(V;.V.-V.V>^ yv'' (A.6) 



82 



since it vanishes identically in an Euclidean space. The Ricci tensor and the curva- 
ture scalar are defined as ^ 



i?"^ -.^ g^"" R'^.p, (A.7) 
V«^;^=[9„52 + ^fJ^;^ y (A.8) 

The metric can also be described through the introduction on every point of 971 of 
an orthonormal basis of d vielbeins {<Jk}k=i- The vielbeins satisfy the properties 



(^kQapCTl = CTakCrt = ^kl 

^1 



aak^ - ^3 (A.9) 



Latin indices are associated to an Euclidean metric. In terms of the vielbeins the 
metric reads 



9al3 = CFakC^pk (A.IO) 

Combined with the compatibility conditions this last equation bares the dependence 
of the Christoffel symbols on the vielbeins. In particular in the Euclidean case when 



dfsaak ^ dacri3k (A.ll) 
the identity holds 

= <Tkd^<^uk = -(J^ikSX (A.12) 
The Christoffel symbol is symmetric in the lower indices by (A. 1 1). 



B Covariant stochastic differential equations 



In general stochastic differential equations are not covariant under change of coor- 
dinates due to the 0{\fdi) increments of the Wiener process. Nevertheless covari- 
ance can be achieved through a path-wise definition of the vielbeins. 
In compact notation the system of stochastic differential equations (25) reads 

^ Here the same definitions of the curvature and Ricci tensor are adopted as in [68,31,5]. In 
[145,159] the curvature tensor has opposite sign but the Ricci tensor as well as the curvature 
scalar are the same as here. 
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g«(T') = Q'" 

g-^iQ') = K4){T') (B.l) 
,(T') = 1 

Geometrically the system describes the transport; of an orthonormal frame of vec- 
tors a^, k — 1, (i parallel to the trajectories of the position process q for any 
given realization of the Wiener process w. In fact, the second equation can be re- 
cast in the form of a covariant derivative along a path q"{t) 



dq'^ = v'^dt+ ^J^a^odwk, 
dat = -^^,(7,^odg^ 
d^ ^-<;fdt, 



IT = ° 

Hence if the metric compatibility condition is imposed at initial time, it will hold 
true all along any given stochastic trajectory. 

In the main text it was stated that solutions of (B.l) are the characteristics of the 
Fokker-Planck equation (22). The statement is verified by differentiating the aver- 
age over the Wiener process 



(<j(T) S^''^ {q{T) - Q)) := J Vfi{w{t)) e't /t' '^^'^(''W'*)^^ (g(T) - Q) (B.3) 

along the trajectories of (B.l). Averages of stochastic differentials are most con- 
veniently performed if the infinitesimal time increment of the Wiener process is 
independent on the current state of the system. Ito stochastic differentials 



a^t) dw\t) lim a," (t) [w\t + dt) - w\t)) (B.4) 

implement the condition [49,91,96,124]. The conversion of Stratonovich differen- 
tial into Ito can always be accomplished by expanding around the pre-point dis- 
cretisation and retaining terms up to order 0{dt) with the proviso 

dwkdwi —Skidt + o{dt) (B.5) 
A straightforward but tedious computation yields 



V 



dat 



hz 

2 m 



2m ^ J \ m 



(B.6) 
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The third equation is discretisation independent. Note that in the Ito equations there 
appear non-covariant quantities. The time differentiation of (B.3) yields 



A(,(T) 5W (g(T) - Q)) = (^(T)<„(g(r),T)-|^5('^) (g(T) - Q) ) 



2m 



-TmnT) 



S^'^ {q{T) - Q)) 



(B.7) 



with the Ito drift 



<to(?m, T) := v^{q{T),T) - ^ r^.(g(T)) g'^''{q{T)) 

Zi Tib 



(B.8) 



The term Unear in dw averages out due to statistical independence of the Wiener 
noise increments. From (B.7) straightforward algebra recovers the Fokker-Planck 
equation. The same result is also obtained by applying functional integrations by 
parts on the original Stratonovich equations [14,122]. 



An analogous calculation evinces the equivalence in measure of the system of 
stochastic differential equations (B.l) with a free Wiener motion on 9Jl 



V m 

da^^-r;^a^,odq (B.9) 
advecting the new potential term 



dg — —Q 



z (j) m I I 

~Y ^ 2hz 



' Z(j) 



1 



WaV" dt + 



v°'aakOdwk 



-g {^dt+\l — v"aak dwk 



(B.IO) 



Equivalence in measure means that the last two equations are associated to the the 
same Fokker-Planck equation as (B.l) 



^/\^\K,(Q,T\Q',T') ^{g(T)5^''^(q{T)-Q)) (B.ll) 

where the average is extended over the solutions of (B.9). The representation (B.ll) 
of the transition probability density is called the Girsanov-Cameron-Martin for- 
mula. 
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Finally it is worth noting that the Euclidean condition (A. 11) is the integrability 
condition for the mapping 



dqk = (^akdq"" 



(B.12) 



which retrieves the natural Euclidean frame where the vielbeins are the versors of 
the Cartesian axes. In the jargon of statistical mechanics, the integrability condi- 
tion permits to map multiplicative into additive noise. In such a case, the Laplace- 
Beltrami operator on a scalar reduces to the Bochner's Laplacian 



g'^'V^V, = [a^, - V^, g''" 9« ] X. = aj^ ( a^d, ) (B.13) 

Hence if (B.12) holds, the Stratonovich equation 



dq'^it) = v''{q{t),t) dt+J— a^{q{t)) o dwk{t) (B.14) 

is covariant since it governs the characteristic curves of the Stratonovich-Bochner 
form of the Fokker-Planck equation [96,97, 165] which is covariant whenever (A. 11) 
holds. 



C Path integrals from stochastic differential equations 



Girsanov's formula (B.l 1) reduces the construction of path integrals to the solution 
of the heat kernel equation on a Riemann manifold Wl for short time intervals. As 
in the main text, 971 is restricted to be either compact and without boundaries or 

W''. The problem can be solved by means of a generalisation of the Levy construc- 
tion generally used to define the Wiener (Brownian) motion in Euclidean spaces 
[104,91]. The derivation which follows is an explicit version of the argument out- 
lined in [50]. 

On R'^, the expectation of the Wiener process at any intermediate time t e [T', T] 
conditioned on its end points is 



The variance is of the order 0{T — T') for t of the order (T — T')/2. The Levy con- 
struction proceeds by dividing the time axis into small slices of duration dT. Within 
each time slice the Wiener process is interpolated with (C. 1). The increments of the 
Wiener process over a time slice are Gaussian random variables 



86 



" dT ^ ^ 

with zero average and unit variance. The projection of the Levy construction over 
a Riemann manifold 9JI gives 



Kit) = -TiMt))<y^{t) (C.3) 

The random equation is now differentiable and equivalent to the geodesic problem 



&) = - r^.(9ge(t)) q',eit) lUt) (C.4) 

The probability distribution of the random system yields the short time expression 
of the heat kernel 

ylg{Q)K,{Q, r + dT\Q', T') - / 5^''\q^,{T' + dT) - Q) (C.5) 

asymptotically in 



0{dT) ~ 0{Q - Q'Y i (C.6) 

The integral on the right hand side of (C.5) requires the solution of the geodesic 
equation (C.4) with the boundary conditions 



g^,(T' + ciT) =g" (C.7) 

On a Riemann manifold the solution is unique if dT is short enough. Straightfor- 
ward algebra gives 



Tfl'^C 2zh 

KJQ, r + dT\Q', T') J (C.8) 



where V is 



Vi, {Q, Q', dT) = aUQ) (^] {T' + dT) 



= - oUQ>0i{Q'){qUT' + c^T), gf,(T')}p.b. (C.9) 
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The Poisson brackets on the right hand side (see appendix E for further details) 
provide the d linearly independent solutions of the linearised dynamics along the 
geodesic [31,68] 



(C.IO) 



Thus det V is the determinant of the scalarised linear dynamics. Since the vielbeins 
are parallel transported along the geodesic 



the scalar fluctuations fulfill the equations of the motion 



(C.U) 



ui{t)=ui{T') it-T') 



T 



One gets into 



(C.12) 



yik{Q,Q',dT) = 6i,kdT 

CTaiiT')a^,{T')R''p,MT')) q^jTXeiT') + o{dT^) (C.13) 

The determinant is computed by means of the identity 



Indet(l-eX) = Trln(l-eX) = -eTrX + o(e) (C.14) 
The approximations hold in probability, thus it is legitimate to replace in (C. 13) 



(C.15) 



which is correct within the leading order in the asymptotics (C.6) [46,99]. The final 
result is 



KMT' + dT\Q'X) 



m 



2 zHtc dT 



2 m\\qge(T')\\^ dT , zhR(Q') 



dT 



(C.16) 
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In order to be elevated to the status of a proof the above derivation requires a precise 
estimate of the errors done in the approximations. The interested reader is referred 
to [1 15,26] for rigorous estimates of the heat kernel on Riemann manifolds and to 
the recent paper [5] for a complete proof of the path integral construction. 

The short time approximation of the propagator is often given in the form [145] 



K,{Q,T'+dT\Q\T')^K,{Q,T'+dT\Q,T')e-fT''^'^o.M.^^^^rp^^ (q17) 
where 



Co.M. ^^\\q-v\\^ + z^(t>+ — Vav'' - ^ — (C.18) 
2 2 12 m 

and [159] 



K,{Q, r + dT\Q, T') = 1 + ^^^^dT + o(dT) (C.19) 

12 m 

In the presence of curvature a rigorous probabilistic result supports the identifi- 
cation of (C.18) as "classical" Lagrangian. Namely for real values of the analytic 
continuation variable z it has been shown [79,151,96,69] that the probability to find 
q{t) in a tube of arbitrarily small radius e around any smooth path r{t) connecting 
Q' to Q in a time interval of arbitrary length [T\ T] reads 



p{\q{t)-r{t)\ < eyt e [T', T]) ~ e-^ir dt^o.M.{r,r) ^^20) 

asymptotically in e. Therefore (C.18) solves the Onsager-Machlup problem of de- 
termining the most probable smooth paths covered by the stochastic process. 
The asymptotics of the diagonal component of the kernel (C.19) describes the 
leading contribution to the square root to the Van-Vleck determinant associated 
to (C.18) 



D Non covariant path integrals 

Beside the covariant formalism discussed above, in the literature are often encoun- 
tered non covariant constructions of path integrals. An example is the Edwards and 
Gulyaev treatment of a quantum mechanical two dimensional particle in a radial 
potential [57,133]. By (B.6) the line element in radial coordinates 
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yields the Ito stochastic differential equations 



, hz , zh , 

ar — — at + \ — awi 
2r \ m 

dO — \ dw2 

\ m r 

zUir) , 

dq ^ ^ dt 

n 



The path integral action in the Ito representation is 



(D.2) 



S^ — / dt— 

zHJt' 2 



T 



m 



r — 



hz 
2mr 



+ r'9 



2 n2 



z h Jt' 



dt 



m 1 



hzr 
~2r 



+ z'U 



(D.3) 



In the Ito representation it is not possible to apply the rules of ordinary calculus. 
In particular, the term in (D.3) linear in the radial velocity cannot be interpreted 
as an exact differential. Only the Stratonovich representation is compatible with 
ordinary calculus. For the term linear in the radial velocity, the relation between the 
two representations is 



hzr 



2r 



hzr 



Ito 



2r 



Strat. 



hz 

imr' 



:{f')dt 



hzr 



2r 



Strat. 



{hzf 

Amr'^ 



(D.4) 



since over infinitesimal increments 



z h 

mdt 



(D.5) 



The first term on the right hand side of (D.4) can be treated as an exact differential. 
The result is the path integral 



x,(i?,e,r|i?',e',r') = J§/ v[r{t)e{tMt)e-^nrr.^^'^' 



m 



r' + r^9^ 



Smr"^ 



+ z'U 



(D.6) 



where the angular kinetic energy is still defined according to the Ito prescription. 
The Edwards and Gulyaev result is finally recovered by setting z — i. 
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Change of variables in path integrals obtained in arbitrary discretisations are no- 
toriously unwieldy. The probabilistic interpretation supplies a useful guideline to 
understand the properties of the path integral. A thorough investigation of the re- 
lation between stochastic differential equations and path integrals can be found in 
ref. [164]. 
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E Summary of some elementary facts in classical mechanics 



The appendix recalls some definitions and results of classical mechanics used in the 
main text. The application of semiclassical methods in quantum mechanics requires 
information not only about the solutions of the classical equations of the motion 
but also about their local and structural stability. While in the first case surveys of 
the material summarised in the present appendix can be found in monographes as 
[7,9,123,110] and [68] for a geometrical point of view, in the author's opinion the 
best introductions to stability problems in classical mechanics remain the classical 
research papers by Gel'fand and Lidskii [71] and Moser [116]. Comprehensive 
presentations of linear Hamiltonian system are available in [163] and [59]. 



E.l Lagrangian versus Hamiltonian classical mechanics 



E.1.1 Lagrangian formulation 

Classical trajectories are extremal curves on some d dimensional manifold 9Jl of 
the action functional 



S = I dtC{q{t)A{t),t) (E.l) 

JT' 

The boundary conditions of the variational problem are determined by some given 
initial conditions. The Lagrangian £ is a function defined on the tangent bundle 
Tm X M of 971. A curve x{t) on TSlJt x R is defined by the Uft of a smooth curve 
q{t) achieved by setting 













dq(t)/dt 



(E.2) 



Velocities transform as well as positions 5" as contra- variant vectors. Lagrangian 
of physical interest are quadratic in the velocities. 

Any smooth extremal of the action must fulfill the Euler-Lagrange equations 



|^_^|^^0 a = (E.3) 
dtdq'' 

The "ellipticity" condition 
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Sp{LjJ > V« 

insures the equivalence of the Euler Lagrange equations to a first order system of 
differential equations in x with fundamental solution specified by a diffeomorphism 
$ [9]. Given an initial condition x' the time evolution of an extremal or classical 
trajectory is then obtained as 



x,e{t) = x',t') (E.5) 

The diffeomorphism ^{t; • , f ) is a flow or one parameter group of transformations 
in t 



M';x',t') ^x' , Vx' 

(E.6) 

<l(t'; x", = ^t; x" , t") , M x" 

The time evolution of infinitesimal perturbations of the initial conditions is gov- 
erned by the linearised dynamics. 



Sx'i.it) = (^1^ j (t ; x', t') 5x'^ := r^it , t' ; x') 5x'^ (E.7) 

The columns of ^ are linearly independent extremals of the second variation along 
a classical trajectory qd'. 



dq°' dt dq'^ 



= {LSqa)„ = (E.8) 



The "classical fluctuations" Sq^ are referred to as Jacobi fields [129,93,51]. The 
linearised dynamics is defined on the space TT^^^QJl tangent to the tangent space 
along the classical trajectory. 



E.1.2 Mother theorem 

Let G be a Lie group with N parameters {t"}^^]^. Acting on configuration space, 
the Lie group generates smooth transformations of variables. In a local neighbor- 
hood of the identity, the transformations are spanned by the N vector fields 



n = l,...,7V (E.9) 

t=o 
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induced by the Lie algebra of G [76,68,120]. The action functional is invariant 
under G if 



r dtC{ip{t),<p{t),t) = r dtC{q{t),q{t),t) 
Jo Jo 

ip{t) = ip{i,q{t)) 
is satisfied. Rephrased in differential form, invariance means 



(E.IO) 



0= 



Vn 



(E.11) 



t=o 



Since Lagrangians of physical interest have the form 



(E.12) 



(E. 1 1) also entails that a symmetry is an isometry of the metric tensor which leaves 
invariant the vector and scalar potential 



a 13 



" dq^ 



+ 9a' 



dqP 



+ 9p' 



" dq-r 
„ dU 



dq^^ 



= 



(E.13) 



Combined with the Euler-Lagrange equations, (E. 1 1) yields 



= 



di 



V 

dq^ " 



(E.14) 



and consequently defines the conserved quantity 



Hr 



dC 
dq^ 



(E.15) 



The index contraction on the right hand side corresponds to a well defined scalar 
product. To wit, under generic change of coordinates q'"' — q'^'iq) the classical 
momentum transforms as a co-vector [7,68]: 



(E.16) 
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E.1.3 Elementary Hamiltonian formulation 

The ellipticity condition (E.4) permits to define the Hamiltonian 



n{q,p,t)^ sup {par-^q,q,t)} (E.17) 

A smooth curve satisfies the Euler-Lagrange equations if and only if it is solution 
of the Hamilton equations in canonical coordinates 



on 

dpa ' 



Pa = 



on 



a — 1, .., d 



(E.18) 



Let $ denote the flow which solves the Hamilton equations. The flow is defined 
in phase space, the geometrical cotangent bundle T*5!Jl. The derivatives of $ with 
respect to the initial positions and momenta specify the flow F of the linearised 
dynamics around a classical trajectory qd- The space where the linearised dynamics 
is defined is tangent to the cotangent bundle TT*^^Tl. A similarity transformation 
connects the linearised dynamics around the same trajectory in phase and tangent 
space: 



f{t,t') = %{t)d{t,t')%-\t') 
by 2 d X 2d dimensional matrix 



(E.19) 



I. 







-qq 



(E.20) 



the blocks being the second derivatives of the Lagrangian evaluated along the clas- 
sical trajectory qd. In particular the lower two blocks specify the pull-back [68] of 
the momentum 



dp{dq) 



Hq)ocf3— + [*-qq)a(30q 



dt 



(E.21) 



In the main text the latter is denoted as 



VSq := dp{6q) (E.22) 
to emphasise that the linear momentum transforms as a covariant vector. 
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E.1.4 Symplectomorphisms 

A general formulation of Hamiltonian dynamics is attained in arbitrary coordinates 
x: 



m 



a , P = 1, ...,2d 



(E.23) 



where f2a/3isa2d x 2d tensor globally defined in T*3Ji and characterised by the 
properties 



— ~ ^/3a , 
detQ 7^ 0, 

dilcep dn.ya dVLjj^ 



+ 



dx^ 



+ 



dx"' 



= 



(E.24) 



In other words p specifies a non degenerate closed two form [3 1,68]. By Poincare's 
lemma [68] any closed two form is locally the curl of a vector potential 



(E.25) 



The elementary formulation of Hamiltonian dynamics is recovered from (E.24), by 
means of Darboux's theorem. This latter states that on a differentiable manifold SDT 
there exist local coordinates x — {q,p) such that the symplectic matrix reduces to 



J := Q 



Daxboux ■ 







-Id 





JtJ = I 



2d 



(E.26) 



A straightforward computation proves the Hamilton equations (E.23) are invariant 
in form under smooth transformations verifying the condition 



(E.27) 



Transformations which satisfy (E.27) are said canonical or symplectomorphisms. 
One parameter symplectomorphisms continuous around the identity are charac- 
terised by the infinitesimal version of (E.27) around the identity: 



*°(a;,t) = x" + V'^{x)i + 0{f) 



(E.28) 
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The condition is satisfied by Hamiltonian vector fields 



= ir^ ^ (E.29) 

the generating function Tiv being a phase space scalar. Canonical transformations 
leave the Poisson brackets of two phase space scalar x and 



invariant. Poisson brackets are antisymmetric, obey the Leibniz chain rule and sat- 
isfy the Jacobi identity [7,110,98] as follows from the properties (E.24) of Q. If 
X and (f) are generating functions the Poisson brackets are equivalent to the skew 
product of the corresponding Hamiltonian vector fields: 



{0,x}p.b. = ^apV^vH (E.31) 

The attribute skew refers to the role of pseudo metric played by the antisymmetric 
tensor Vt^js. The skew product vanishes on pairs of linearly dependent vectors. As a 
result in a phase space of 2 d dimensions the maximal number of linearly indepen- 
dent vectors which can be pairwise annihilated by the skew product is d [7,1 10]. 



E.1.5 Symmetries 

Poisson brackets supply a convenient formalism to handle the relation between 
continuous symmetries and conservation laws. Namely if 7^^, TYj, are respectively 
the generating functions of the one-parameters symplectomorphisms ^a(") t") and 
^^(•) t'')' the derivatives along the flow obey the chain of identities: 



' &Ha 



(*6(x,t')) 



t''=0 



(E.32) 

t"=0 



An integral of the motion is therefore the generating function of a symplectomor- 
phism leaving the Hamiltonian invariant. In general, Frobenius theorem establishes 
integrability conditions in terms of the Lie algebras of vector fields [68]. Poisson 
brackets are an integral version of the Lie brackets of the Hamiltonian vector fields 
with generating functions Ha,Hb [7,110]: 



{Va , V,}U_ := - V,^^ = - ^^"^^{^a , ^Jp.b. (E.33) 
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If the generating functions Poisson commute mutually and with the Hamiltonian 
7i they are simultaneously preserved by the dynamics. The corresponding Hamil- 
tonian vector fields form an Abelian group of transformations. More generally, a 
symmetry group is realised by a Lie algebra of Hamiltonian vector field {K}^i 
with structure constants C^^ (see appendix F) 

{K,H}£.b. = Cl,V^ (E.34) 
yielding the Poisson bracket algebra 

{Tia , Tibjp.b. = -Cl^Hc + Dah, Dab = - (E.35) 

The -Dab's are skew symmetric phase space constant depending only on the struc- 
ture of the Lie algebra [7]. More light on the meaning of the Dah^ is shed by 
the analysis of the dependence of the {Ha}a=i the Lie algebra vector fields. 
In consequence of (E.24), (E.25) the generating functions must admit the general 
representation 

^a-^aC + ^a (E.36) 

for {ha}a=i some scalar functions and the vector potential defined by (E.25). 
The definition of Lie brackets (E.35) entails the identity 



UVa , = 1^ - 1^ - ^^}p.b. (E.37) 

which, combined with (E.34) and (E.36), yields 



Da.-Cl,K-Va^^^+V,^^^ (E.38) 

Thus it is only when the -Dafc's are zero for all a, h that the scalar functions {ha}a=i 
can be set to zero in the representation of the generating functions of the Lie alge- 
bra. In such a case a symmetry is said equivariant. The denomination is justified 
by considering an example with globally defined canonical coordinates. In such a 
case the vector potential is 



= (E.39) 
while Hamiltonian vector fields have the form 



Va - (E.40) 
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Hence one has 



^aC = ^« + - 2 Ha^ (E.41) 

Equivariance then means that the constant of the motion is homogeneous of degree 
two in X. The rotation group provides an example of equivariant symmetry leading 
to the conservation of the angular momentum 



Tia^Pmrl^iqu (E.42) 
with r(^)^ the infinitesimal rotation matrix around the a axis. 

E.2 Linear Hamiltonian systems 
E.2.1 General properties 

Poisson brackets allow a covariant description of the classical linearised dynamics 
in the space TT*^^9Jt tangent to phase space along a classical trajectory 

Xct{t)^m]x',t') (E.43) 
The linear evolution matrix can be recast in the form 

v-,{t^t')^n^\x')n,,{x') ^ (^^J {t,x\t') 

= , <,(t)}p.b. ^p^{x,i{t')) (E.44) 

In matrix notation, the linearised flow satisfy the linear Hamiltonian equations 

^{t\t') = hd (E.45) 
having defined 



dx'^dxP dx^ dxP 



(E.46) 
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(E.44) provides for the invariance of the skew product of any particular solutions 
of the linearised dynamics under canonical transformations. Namely any canonical 
transformation of the flow 



a^ = ^(-^^W)F^(M') (E.47) 
satisfies the symplectic property 



Flit, t') n^sit) Fj(i, t') = n^fsit') (E.48) 

Combined with (E.44), the equality specifies the behaviour of the linearised dy- 
namics under time reversal 



if-'rpit,t') = F^(t',t) (E.49) 

The property is inherited by the linear evolution matrix ^ in TTg^^Tl. 

By (E.32) symmetry transformations ^ commute with the Hamiltonian flow $ 



o ; x', t') = ^{x'),t') (E.50) 

Hence the linearised dynamics must satisfy 



F^(t,i';^(a:')) = ^MWm';x')^{x') (E.51) 



E.2.2 Linear periodic systems 

Linearisation around a classical trajectory on a periodic orbit of period Tee gives 
rise to a periodic Hamiltonian matrix 



H{t) = H{t + Ta), yt (E.52) 

The pseudo metric has the same periodicity. The general form of the solution 
of a linear periodic system is dictated by Floquet theorem [9,163]: 



F(t,t') = Pe(t , t') exp <^ —- / ds (l^"^H)(s) \ (E.53) 



Tct Jt' 

with Pe a periodic matrix such that 
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Pe{t' + nT,i;t') ^ \2d (E.54) 
for all integer n. 

The stability of a classical periodic orbit is governed by the monodromy matrix: 



M{Ta) := F{t' + T,e,f) = exp 



{ir 



ds (Jl"^H)(s) 



(E.55) 



It is not restrictive to choose local coordinates such that the monodromy matrix 
satisfies 



Mt JM = J 



Together with 



detM = 1 



(E.56) 



(E.57) 



(E.56) defines the linear symplectic group 5*^(2 d). Any symplectic matrix has the 
square block form 



M 



A 
C 



B 
D 



Dt 

-ct 



-Bt 
At 



(E.58) 



Furthermore the symplectic condition (E.56) requires the d x d blocks to fulfill 



ADt - BCt = \a 

ABt = BAt , CDt = DCt 



(E.59) 



or equivalently 



DtA-CtB = Irf 
DtB=BtD, CtA = AtC 



(E.60) 



Every linear Hamiltonian flow in Darboux coordinates draws a curve in the sym- 
plectic group. The normal forms of elements of Sp{2d) are also strongly con- 
strained by (E.56). Left (generalised) eigenvectors of a symplectic matrix are spec- 
ified by the right (generalised) eigenvectors through a complete set of skew orthog- 
onality relations [116,35,59]. In this way it is possible to construct a symplectic 
basis with elements {Zn, U}t=i satisfying 
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-'m,n 



(E.61) 



which reduces simultaneously the symplectic matrix and the pseudo metric J to 
normal form. 

The eigenvalues of any element of Sp{2 d) are constrained in a rigid pattern of pairs 

or quartets in the complex plane. As a consequence the normal form of a symplec- 
tic matrix typically (but not necessarily!) consists of the two or four dimensional 
blocks [116,35,163] listed below 

1 Direct hyperbolic blocks 
The eigenvalues are (e'^ , e~'^) with uj real. The simultaneous normal forms of 
an hyperbolic block and of the associated pseudo-metric are: 



Md.i 











Jd.h. — 



-1 



(E.62) 



2 Inverse hyperbolic blocks 

The eigenvalues are ( — e'^, , — e~'^), a) real and 



Mi, 



- 









Ji.l 



-1 





(E.63) 



Inverse hyperbolic blocks occur because the monodromy matrix may not have 
a real logarithm for all odd iterates [100]. 
3 Elliptic blocks 

The eigenvalues are (e"^ , e~ *'^) with u real. The normal forms are 



Mell. = R2(-C^) , 



cos a; 
- sin a; 



■Jell. - 

sin a; 
cos a; 



dR2 
djjj 



(0) 



The eigenvectors 



(E.64) 



(E.65) 



are at the same time orthogonal with respect to the skew and the standard 
scalar products. Thus they satisfy 



sign 



(E.66) 
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The signature of the skew products is a constant of the motion known as Krein 
invariant. It provides an intrinsic characterisation of "positive" and "negative" 
frequencies [95,71,116,59]. 
4 Loxodromic blocks 



The eigenvalues are (e'^+"^ , e' 



M 



lox. 



Jlox. 



e^R2{-cu) 

e-^R2(-u;) 

-I2' 



I, 







(E.67) 



In a one parameter family of real linear Hamiltonian systems the creation of 
a loxodromic block requires the simultaneous exit of four eigenvalues from 
the real or imaginary axis. Such bifurcations will occur only at isolated values 
of the parameter corresponding to points where two pairs of hyperbolic or 
elliptic eigenvalues coalesce. 
Parabolic blocks 

A parabolic block comprises two degenerate unit eigenvalues paired with an 
eigenvector and a generalised eigenvector: 



M 



pax. 



■"par. 







(E.68) 



The non diagonal entry k is fixed by the skew product of the elements of the 

basis yielding the normal form. 

From the differential equations point of view, a parabolic block is the value 
over (a multiple of) the period of a marginally unstable block of the linear flow 
F. Such block is spanned by linear combinations of periodic vector fields with 
coefficients polynomial in time [129,51]. A simple eigenvalue degeneration 
corresponds to the evolution in phase space of a periodic Xi (t) and a non 
periodic X2 (t) Jacobi fields 



f{tX)x,{r)^x,{t) 

F{t, T') X2{T') = X2{t) = K l^^xi{t) + y{t) 



The conservation of the skew product 

\ = {x\^xx){t), Ml 



(E.69) 



(E.70) 



enforces the normalisation condition whence (E.68) follows. Parabolic blocks 
are a generic feature of systems with continuous symmetries. Jordan blocks 
also appear in the presence of unstable degenerate eigenvalues different from 
unity [129,71,116]. 
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Hyperbolic and loxodromic blocks characterise unstable directions of periodic or- 
bits. Parabolic blocks are marginally unstable and exhibit a linear growth of the 
perturbation along the direction spanned by the generalised eigenvector. Finally 
elliptic blocks describe under certain conditions stability [71,1 16]. If the character- 
istic frequencies {un}n^'^ of the elliptic blocks are mutually irrational the linear 
flow is stable. Furthermore, Moser has proven in ref. [116] that a periodic orbit is 
almost stable versus nonlinear perturbations if the linearised flow is elliptic with 
mutually irrational frequencies. Almost stable here means that there exists a formal 
power series (possibly divergent) supplying a Ljapunov function for the periodic 
orbit. In consequence a parametric perturbation of the periodic orbit will remain 
for extremely long times in the neighborhood of the periodic orbit. 

Finally it is worth stressing that a block contributing to the monodromy with a 
certain stability may change stability at generic times [59]. As a matter of fact, the 
symplectic property requires the eigenvalues of a linear Hamiltonian flow to be con- 
tinuous under time evolution but does not rule out discontinuities of the derivatives. 
Inverse hyperbolic blocks are generated by the time evolution from bifurcations of 
unstable elliptic blocks. The phenomenon is exemplified by the dynamics of a lin- 
ear Hamiltonian system with Hamilton matrix periodic and positive definite. Rule 
out the trivial example of the harmonic oscillators and assume: 



- [Sxl^ H Sxci] = Sxi^ — Sxct > , (E.71) 

At initial time the eigenvalues are equal to one. Since the assumption (E.71) forbids 
the formation of hyperbolic blocks at unity, as time increases the eigenvalues can 
only move on the unit circle. In particular Krein positive eigenvalues (E.66) move 
counterclockwise on the upper half unit-circle and the Krein negative ones on the 
lower half-circle. If two eigenvalues with opposite Krein signature meet at minus 
unity, they can satisfy the symplectic condition in two ways. They can cross each 
other and continue their motion on the unit circle or they can leave the unit circle 
and start moving on the negative semi-axis. Note that the second option is ruled out 
for eigenvalues with the same Krein signature. A bifurcation to inverse hyperbolic 
would imply in that case a change of the overall Krein signature. The eigenvalues 
will at some later time come back to minus unity and resume their motion on the 
unit circle, with Krein-negative ones now moving clockwise on the upper-half cir- 
cle. All of that can happen for a stable T^^-periodic system: the eigenvalues must 
only be back on the unit circle at times multiple of the prime period. More infor- 
mation can be found in [71,116,59,163]. 



E.2.3 Eigenvalue flow around a parabolic block 

Consider a one parameter r family of linear Tc^-periodic Hamiltonian systems 
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(E.72) 



with analytic dependence in r. In what follows at r equal zero the parametric de- 
pendence will be simply omitted. At fixed times the linear flow evolves in r ac- 
cording to a linear Hamiltonian equation 



£dt'F;'^{t,t')6U{t') f;\t,t') frit,T) (E.73) 



with 



SH{t) = Hr{t) - H{t) (E.74) 

The parametric stability is determined by the r dependence of the eigenvalues of 
the monodromy matrix. Up to first order accuracy in r, this latter is 



= M + r ^^j^^^^dt F-^\T+Ta, t) 
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f-\T+T,e,t)M (E.75) 

r=0 



with 



J = Q(T) = n{T + Tee) (E.76) 

The leading correction to the eigenvalues is obtained projecting the above equation 
on the unperturbed eigenvectors. 



An isolated zero mode of the periodic Sturm-Liouville fluctuation operator pairs 
up with an elementary parabolic block of the monodromy as (E.68). Multiple zero 
modes brought about by Abelian symmetries also decouple to form an equal num- 
ber of elementary parabolic blocks. For generic 5H, the qualitative effect of the 
perturbation in r in such cases is to shift away from unity the eigenvalue pair of 
the parabolic block. Since the symplectic structure is preserved, a parabolic block 
generically bifurcates into an elliptic or hyperbolic one. Groups of unit eigenvalues 
of the monodromy corresponding to non- Abelian symmetries or marginal degener- 
ations can also be analysed in a similar way [71]. 

The right periodic eigenvector is the phase space lift of a periodic Jacobi field, 
a zero mode of Lpcr. in [T,T + T^^^]. If the zero mode stems from a continuous 
symmetry according to Nother theorem in Darboux coordinates it will have the 
form 
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V(5t(Zc£(T) 



(E.77) 



SiQce is the vector field induced by the infinitesimal generator of the symmetry 
transformation. Together with a generalised eigenvector t2, (E.77) specifies a dual 
basis: 



Mti = ti 

M r2 = r2 + K Ci 

whence it follows 



Mt J n = J ti 

J r2 = J r2 - K J ti 



(E.78) 



provided 



[[ti = 4Jri 



(E.79) 



(E.80) 



Here, the class of perturbations of interest comprises those ones which do not 
change the Morse index of the Sturm-Liouville operator associated to the second 
variation of the Lagrangian around the periodic orbit. Any strictly positive definite 
Tc^-periodic term \J{t) added to the potential in the Sturm-Liouville operator 



6q'< Lqq-rSq = 6q^[lqq + T \J] 6q 
implements the condition: 



(E.81) 



dt{dq^\Jdq){t) > 0, 



W6q{t) 



(E.82) 



The perturbation of the Lagrangian introduces in Darboux {q , p) coordinates the 
Hamiltonian perturbation 



5H{t;T) = T 







(E.83) 



Up to leading order the eigenvalue equation for the perturbed parabolic block re- 
duces to 
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= det 



1 — m(r) K 
tV l-m(r) 

with V specified by 



(E.84) 



V: 



dH 



f-\T+T,e,t)M 

T=0 



(E.85) 



The eigenvalues are 



m±(T) = 1 ± K T/| e^'^ ^ + 0(t) 

~ exp <^ ± K 1^1 e'"" ^ | + 0(r) 



(E.86) 



The exponentiation is legitimate because of the constraint imposed by the symplec- 
tic structure. It determines the phase factor in the exponential modulo 2 tt. 

There are two main lessons to be drawn from (E.86). First, the characteristic fre- 
quency has a power expansion in ^/T and not in r as the monodromy matrix. This 
is not surprising as the frequency of an harmonic oscillator appears quadratically in 
the equations of the motion. The eigenvalue split induced by an Hamiltonian per- 
turbation on blocks of higher eigenvalue degeneration N is amenable to the power 
series expansion [71,18] 



k = 1,...,N 



(E.87) 



The second lesson is that the nature of the block emerging from the bifurcation 
is invariantly determined by the signature of the product V k. According to the 
conventions adopted for the symplectic pseudo-metric, V is positive definite and 
the signature of the non-diagonal element in the parabolic block directly specify 
the bifurcation. 

The physical meaning of the non diagonal element in the parabolic block is better 
realised by considering the "default" zero mode encountered in the Gutzwiller trace 
formula. 

By energy conservation classical periodic orbits occur on closed curves of level 
of the Hamiltonian E = 7i(x). Periodic orbits are then expected to appear in one 
parameter families the period Tci{E) whereof smoothly depends on the energy. In 
such a case the periodic eigenvector of the monodromy matrix is the phase space 
"velocity" of the classical trajectory 



(E.88) 
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In the absence of further conservation laws Xce{t) is generically the only periodic 
eigenvector in [T,T + Td]. The associated generalised eigenvector is found by 
observing that for all t: 



Xci{t + Ta,{E, ...) ; E, ...) = Xce{t ; E, ...) (E.89) 

The derivative with respect to the energy yields the phase space lift of a non periodic 
Jacobi field. One can identify in (E.69) 



(E.90) 

-rp — =constant 

the second vector field being periodic by construction. Thus y{T) provides the 
sought generalised eigenvector X2. Furthermore energy conservation along the tra- 
jectory 



u\ ^^^^ (+\ 

^2W = ^W> 



dx 



dE 



it) 



^ ^ (m) ^^^^^^^^ ^ (^'^*^^) ' ^^-^^^ 
yields for the family of periodic orbits 



KiE) = -^{E) (E.92) 

Hence the non diagonal element of the parabolic block is equal to the variation in 
period of two orbits, infinitesimally separated in energy. 
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F Lie's first fundamental tlieorem 



A Lie group is a smooth manifold G on which the group operations of 
product : G xG g(t) ■ g(s) = g(/(t,5)) 

(F.l) 

inverse : G ^ G g(t) • g(t(-^)) = g(0) = identity 

are defined. The analytic mapping / = (/^, .., f^) governs the composition law 
governing the group operations 

t" = r(t,o) = r(o,t) 

r(t,/(s,t))) = r(/(t,5),o) (F.2) 

A Lie group acts on a d-dimensional configuration space 971 through a smooth 
mapping (p 

: G X SOT ^ mi (F.3) 
The mapping induces transformation laws of point coordinates of 2Jl 

?[^] = <^"(?,t) (F.4) 

The origin in the {tj-jj_^ space is chosen to correspond to the identity transformation 
onSW 

99"(g,0) = g« (F.5) 

Lie's first fundamental theorem [76] relates the derivatives of (F.4) at a generic 
point t of G to the vector fields induced on 9JI by the infinitesimal generators of the 
group transformations. 



- (F6) 

t=o 



The group composition law (F.2) permits to use both left or right infinitesimal trans- 
lations at t. A left translation at t is 



(^"((^(g,t),5) = (^"(g,/(5,t)) (F.7) 
Differentiating both sides with respect to the s's in zero yields 
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fi=0 



(F.8) 



The matrix 



(F.9) 



s=0 



characterises the infinitesimal left-translation. 

Analogously, exchanging t with s in (F.7) it is possible to express the derivatives of 
group transformations at t in terms of an infinitesimal right translation 



s=0 



^(,,t):=(9l-)^(t)^(,,t) (FIO) 



Comparing the two expressions (F.8), (F.IO) one finds 



%{q, t) = t) (g) = £^(t) t)) 



at' 

and therefore 



(Fll) 



(iH-^)^(t)£g(t)<((^) 



(F.12) 



This latter equality specifies the adjoint representation of the action of the group. 
Namely a right translation can be represented in the guise of a left translation: 



ip'^iipiq, B),t) = ip'^iipiq, t), fit, /(S, t(-^))) 
exploiting the expression of the identity 



(F13) 



Differentiation at s equal zero yields: 



(F.14) 



A<(t)=^(t,/(B,t(-)) 



s=0 



(F.15) 



From the representation (F.ll) of the derivatives of group transformations it is 
straightforward to derive the structure constants of the group. To wit, the existence 
of a global parametrisation of group transformations in terms of the variable t re- 
quires 
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The identity implies for infinitesimal left translations 



(F.16) 



2," £m 



-a' 13 



dvf 




'82," 


821] 








9t« 



(F.17) 



The vector fields do not depend explicitly on the group coordinates t. Hence the 
condition (F. 17) admits solution if and only if it is possible to separate the variables. 
In other words there must exist some constants ^ such that 



/^c ria' rife' 



(F18) 



is satisfied. The constants C^^ are the structure constant of the group while (F.18) 
are the Maurer-Cartan structure equations. 

Finally, in order to prove that (292) is the invariant measure of the group one ob- 
serves that if 



VP"(g,t) = v^"(<^(g,0),5) 
differentiating the first equation with respect to yields 

Therefore one gets into 



(F.19) 



(F20) 



or equivalently 



(F.21) 



whence it follows that 



(F22) 



N 



dG=l[ de det £(t) 



(F23) 



a=l 
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is the right invariant measure. However for a compact connected group the invariant 
measure is unique (up to a constant factor) and therefore also 

det£(t) ocdet9^(t) (F.24) 
must hold. 

More details can be found in [68,76,120,144]. 
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G Fresnel Integrals 



The paradigm of Fresnel integrals is provided by the one dimensional case 



^l-^e'^ (G.l) 

with z a real number. The integral is not absolutely convergent since the integrand 
is in modulo equal to one. Nevertheless intuitively one can hope that the integral 
converges on the real axis due to the increasingly fast oscillations of the integrand. 
A quantitative analysis can be performed on the complex g-plane. Since the in- 
tegrand is even, it is enough to consider the first quadrant of the complex plane, 
R+ X i R+. The integral can be made absolutely convergent if 

Re^izq^} < ^ cos (2 arg + arg 2; + < (G.2) 
A Gaussian integral is recovered each time 



TT 

2 arg 5 + arg 2; + - = TT 

The absence of poles in the domain of convergence in the complex plane permits 
to enclose the Fresnel integral into a null circuit receiving the other non-vanishing 
contribution from a Gaussian integral: 



dqe 2 = / dqe"^ 2 — e*4 / d\q\e 2 if argz = 
Jk+ J'r+ 

^ j) dqe'^ ^ e~'^ J d\q\e^ - j dqe'^ if arg^; = 7r(G.3) 



The final result is 



dq jzq^ 



dz) = / . — 
The multidimensional generalization is 



27r 

6 4 



(G.4) 



d^q lW, e-f-'i-'-^"' 

nre' ^ = , (G.5) 



ind~L^^) being the number of negative eigenvalues of the symmetric matrix L^^^ 
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Quadratic path integrals are the continuum limit of a lattice Fresnel integral ob- 
tained from the discretisation of the action 



n 

defined by mid-point rule 

4^^ = \ {mn) - Sq{n - l)]t L^^(n) [Sq{n) - Sq{n - 1)] 

^rr / \ r / \nt , / . Sq(n) + Sq(n — l) 
+2 [6q{n) - 6q{n - 1)]^ L<j,(n) ^ ^ 

, [(5g(n)+^g(n-l)]t , ^ + - 1) | 

2 2 



(G.6) 



(G.7) 



with 



5g(n) = Sq{T' + nAt) , 



T-T' 

N 



The continuum limit does not depend on the discretisation of the potential term. The 
form (G.5) of the Fresnel integral is attained by collecting the quantum fluctuation 
in a single A^^d-dimensional vector AQ: 



= AgtLW AQ 
SQ^ = \5q\0) , dq'iO), ...,6q\N),..., Sq^N) 



(G.8) 



The matrix L^^^ depends both on the discretisation and the lattice boundary condi- 
tions. For example, a one dimensional harmonic oscillator with periodic boundary 
conditions yields 



2 + -1 
-1 2 + -1 



-10 ... -1 2 + uj^ 
From (G.6) it follows that the lattice path integral is finally 



(G.9) 



AN) 



d'^5q{n) 



(2 7rnAi)<^/2 e*<^i 



(G.IO) 
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H Some exact path integral formulae 



The stationary phase is exact for quadratic integrals. Namely its effect is to decouple 
the classical from "quantum" contributions to the action functional. The propagator 
path integral becomes 



K{Q,T\Q',T') = / V[Sq{t)]e 



S{Vh&q) 

h 



^e'^L{T,T') (H.l) 

All the spatial dependence is stored in the action function evaluated on the clas- 
sical trajectory matching the boundary conditions. The path integral l{T,T') 
reduces to a pure function of the time interval. The observation allows to shortcut 
the analysis of the continuum limit by a self consistency argument. 



H.0.4 Free particle propagator 
The free propagator is known to be 



p 2h (T-T') 4 

K,UQ,T\Q',T') = (H.2) 

(2 TT n.) 2 

Comparison with (H.l) fixes the normalisation of the path integral to 



Kfree(0,r|0,r') = tfree(T,r') (H.3) 

In consequence the general formula becomes 



K{Q,T\Q\T') = e^%^ Xfree(0, r|0, T') ^^^^ (H.4) 



iiree{T,T') 



H.0.5 One dimensional harmonic oscillator 
The evaluation of the action functional 



s= r dt 

Jt' 



2^ 



(H.5) 



on a classical trajectory connecting Q' to Q in AT = T — T' (open extremals) 
yields 
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Sc,{Q,T\Q',T') = 



mu! 



2 sin(c; AT) 
Divergences are encountered at 



iQ^ + Q'^) cos(a; AT) -2QQ' (H.6) 



ujAT = nn 



(H.7) 



The divergences signal the existence of periodic orbits. The quantum propagator in 
such cases reduces to a Dirac S function. 

Quantum fluctuations around open extremals are governed by the Sturm-Liouville 
operator 



5q{T') = 5q{T) = 



(H.8) 



The operator is diagonal in the complete, for Dirichlet boundary conditions, basis 
of odd Fourier harmonics 



Xn{t) 



AT 



sm 



nt 
AT 



n 



1,2, 



(H.9) 



On a time lattice an infinite dimensional orthogonal matrix changes the variables 
of integration in (H.l) from the quantum fluctuation Sq to the amplitudes of the 
Fourier decomposition. The effect of the quantum fluctuations is enclosed in the 
eigenvalue ratio 



1 



c^AT^ 



rnr 



i^ind-Lui,.{[T',T]) 



e 2 



sin (a; AT) 



CO AT 

The last equality follows from the analytical continuation of the ^-function 



(H.10) 



as) = sr=o 



Res > 



(H.11) 



The Morse index is the number of negative eigenvalues of the Sturm-Liouville op- 
erator (H.8) 



ind-LDir.([r', T]) = # ^ n I 1 - ('^^^ < ^ = Int 



n 



V) 



a; AT' 



TT 



(H.12) 
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The trace of the quantum harmonic oscillator propagator corresponds to the integral 



Tr K{T,T') 
^ [dQ 



2 7rh |sin(u;AT)| 



1 



2 sin 



a; AT 



g-^ [ "^°""'^^^^+V ind-LDi,.([r',r])] 



(H.13) 



The trace is also the inverse square root of the determinant of the self-adjoint oper- 
ator 



■ UJ 



6q{T') = 5q{T) , 



dt ^ ' dt ^ ' 



(H.14) 



The latter admits as eigenfunctions both odd and even harmonics of the Fourier 
basis in the interval [T', T]. The path integral measure becomes 



V\5a(t)\ -^S^n ^C-^-^^c^'^e-^ 
Thus the eigenvalue spectrum of (H.14) 



(H.15) 



27rn 



n = 0,±l ±2,... 



yields immediately the Morse index 



(H.16) 



ind-LperXfT'', T]) = 1 + 2 Int 



27r 

sign sin (u; AT) 1 



+ Int 



cuAT 



TT 



(H.17) 



while the absolute value of the determinant can be extracted as above from the 
C-function. 
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